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ABSTRACT 

Surrounding SN 1987 A is a three-ring nebula attributed to interacting stellar winds, yet no model has success- 
fully reproduced this system. Fortunately, the progenitor's mass-loss history can be reconstructed using light 
echoes, in which scattered light from the supernova traces the three-dimensional morphology of its circumstel- 
lar dust. In this paper, we construct and analyze the most complete map to date of the progenitor's circumstellar 
environment, using ground and space-based imaging from the past 16 years. PSF-matched difference-imaging 
analyses of data from 1988 through 1997 reveal material between 1 and 28 ly from the SN. Previously-known 
structures, such as an inner hourglass, Napoleon's Hat, and a contact discontinuity, are probed in greater spa- 
tial detail than before. Previously-unknown features are also discovered, such as a southern counterpart to 
Napoleon's Hat. Careful analyses of these echoes allows the reconstruction of the probable circumstellar en- 
vironment, revealing a richly-structured bipolar nebula. An outer, double-lobed "Peanut," which is believed 
to be the contact discontinuity between red supergiant and main sequence winds, is a prolate shell extending 
28 ly along the poles and 1 1 ly near the equator. Napoleon's Hat, previously believed to be an independent 
structure, is the waist of this Peanut, which is pinched to a radius of 6 ly. Interior to this is a cylindrical hour- 
glass, 1 ly in radius and 4 ly long, which connects to the Peanut by a thick equatorial disk. The nebulae are 
inclined 41° south and 8° east of the line of sight, slightly elliptical in cross section, and marginally offset 
west of the SN. From the hourglass to the large, bipolar lobes, echo fluxes suggest that the gas density drops 
from 1-3 cm -3 to > 0.03 cm" 3 , while the maximum dust-grain size increases from ~ 0.2/im to 2fim, and the 
silicatexarbonaceous dust ratio decreases. The nebulae have a total mass of ~ 1.7M Q . The geometry of the 
three rings is studied, suggesting the northern and southern rings are located 1.3 and 1.0 ly from the SN, while 
the equatorial ring is elliptical (b/a < 0.98), and spatially offset in the same direction as the hourglass. 

Subject headings: circumstellar matter — dust — scattering — stars: mass loss — supernovae:individual (SN 
1987A) — techniques: image processing 



1. INTRODUCTION 

Supernova (SN) 1987 A in the Large Magellanic Cloud 
(LMC) is the nearest SN in 400 years; the only SN for which 
the progenitor star has been identified and studied; and the 
only SN whose evolution is resolved and observed at all wave- 
bands. Located approximately 50 kpc away, this system is just 
marginally resolved by most ground-based telescopes, even 
with adaptive optics and high-resolution cameras. Matters 
are further complicated by bright, neighboring stars, which 
severely crowd the innermost field of interest for ground- (and 
even some space-) based imaging and spectroscopy. Nature 
has a terrible sense of irony in placing the first nearby SN 
since the inventions of the modern telescope and detector at 
the far extreme of their observing limits. 

The progenito r has been identifie d as Sk -69° 20 2 
(Sanduleak ll969l) . a B3 I supergiant ( Walbornet alJ il989). 
From its spectral type and distance, one infers a luminosity 
of ~ 10 5 L Q , a surface temperature 15000-18000 K, and a 
radius around 43/?©, whic h imply that the star exploded as 
a blue supergiant (or BSG: fWooslev et aflll987l) . This small 
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radi us is co nfirmed by the underluminous light curve of the 
SNlA rnettl 19871) and the extremely-short time delay between 
dete ction of the neutrino burst and o ptical brightening of the 
SN IShigevam aet alll987tlArnettll988l) . 

Sk -69° 202, has a highly-structured circumstellar environ- 
ment (CSE), evidenced in large part by the famous three-ring 
nebula (shown in Fig.Q. This material was flash-ionized by 
the SN light pulse, and is observed today through recombina- 
tion cooling. The kinematics of the inner, equatorial ring (ER) 
show it to be a planar ring, expand ing at ~ 1 . 3 km s~ rathe r 
than a limb-brightened ellipsoid (Crotts & Heathcote 1991). 
Assuming the ER is c ircular, its observed geometry i mplies 
it is inclined at -43° JPlait et alJll995tlBurrows et alJI 19951) . 
with the northern edge closest to earth. The morphology of 
the outer rings (ORs) is poorly understood, in part because 
the south OR (SOR) is distorted from the elliptical shapes of 
the ER and north OR (NOR). 

While SN 1987A exploded as a BSG, the presence of 
circumstellar material, as well as significant CNO process- 
ing within that gas, indica tes that the star first passed 
through a RSG phase tFransson et alJ 119891) . To form 
the three rings, many authors have looked to interact- 
ing stellar winds (IS W) between the RSG and BSG out- 
flows . In this model llKwokl Il982t iBalick. Preston & Ickd 
1987), an equatorial overdensity in a previously-expelled 
slow, dense RSG wind focuses a fast, tenuous BSG wind 
into a polar trajectory, the interaction from which pro- 
duces a bipolar peanut-sha ped neb ula, or wind-blown bubble 
JWooslevll988llWooslev. Pinto & WeayeJll988tlArnett et alJ 
1 19891 IChevalier & Emmerind 119891 iFransson et all 119891: 
ILuo & McCravl1991hH WanP & Mazzal lT99l . Here, the ER 
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FIG. 1.— ffiTWide Field and Planetary Camera 2 (WFPC2) negative im- 
age of a 4',' 5 X 5"0 field surrounding SN 1987A, taken in F656N. North is 
up, east is left. The ER (central ring) surrounds the ejecta from the SN (cen- 
ter), and is flanked by the north and south ORs (NOR and SOR, respectively). 
The co mpanion stars 2 and 3, using the classification of Walker & Suntzefi 
1 1990), are positionally coincident along the line of sight. To increase the 
display range, the ER is displayed with a separate color stretch. 



is the overdense waist of the bipolar nebula, and the ORs lie 
at some point along its outer edges. Such a structure of dust 
has been traced out by scat t ered-li ght echoes as reported in 
ICrotts. Kunkel & Heathcotel i 19951 hereafter CKH95), sup- 
p orting this general mod el. 

iBlondin & Lundavistl (fl993) first modeled the ring nebu- 
losity within a full 2-D, time-dependent hydrodynamic re- 
alization of the ISW model, assuming reasonable values for 
the mass-loss rates and expansion velocities of the stellar 
winds in both RSG and BSG stages, as w ell as an ad hoc 
asymm etry function for the RSG wind. Martin & Arnett 
( 1995) also performed this exercise with a simple cooling 
law, creating synthetic images from their model (using an 
assumed opacity and ionizing UV flux) which reproduce 
the ring geometry only if one identifies the ORs with limb- 
brightening along peanut-shaped bubbles rather than discrete 
rings. Rather than rely on ad-hoc asymmetry functions, 
ICollins et alJ i ll 9991) incorporated the wind-comp ressed disk 
theory (or WCD, Biorkman & Cassinelli 1993), in which 
winds leaving a rapidly rotating star converge at the equator 
as they orbit the star, to provide the equatorial RSG over- 
density. They find good agreement with the ER geometry 
but require that the star rotate at 30% of the breakup rota- 
tion rate. To spin the envelope this quickly, they invoke a 
binary merger in which a companion spun up the envelope, 
then joined the primary core llLivio & S oker 1988). None of 
these models have provided a satisfactory explanation of the 
ORs, leading some to conside r in more detail the influence of 
a binary companion ( Podsiadlowski, Fabian, & Stevens 1991; 
iLlovd. O'Brien. & Kahnll995l) . 

These are only a few examples of a much larger body of 
work, and despite these extensive efforts, many of the critical 
details of the final evolution of the progenitor and of the CSE 
surrounding the SN remain speculative. Within the generally- 
accepted framework of the ISW formation hypothesis lie a 
large number of variant models, invoking magnetohydrody- 
namics, binary companions, multiple evolutionary tracks, di- 
verse mass-loss histories and dichotomous final states. Un- 
fortunately, no m odel has accurately reproduced this system 
( Sugerma neit all2005l) . 



Thus far, formation models have assumed as sufficient con- 
ditions a set of "reasonable" (i.e., typical of nearby LMC su- 
pergiants) values for the mass, radius, temperature and com- 
position of Sk -69° 202, M and v ap for the RSG (1 -2 x 10" 5 
Mq yr\ 5-20 km s" 1 ) and BSG (10 -6 - 10" 7 M yr" 1 , 300- 
600 km s -1 ) winds, and have used as the primary constraints 
the few physical conditi ons of the ER of w hich we are cer- 
tain: the radius (0"81; Burrows et al.lll99 5[). and th e expan- 
sion velocity (10.3 km s" 1 ; Crotts & Heathco tel 1 99 ll) . Results 
matching only a handful of observed conditions can only be 
regarded as suggestive and preliminary. However, up to this 
point, only a handful of constraints exist. Our understanding 
of the circumstellar history of SN 1987 A is limited by a lack 
of observed constraints, not theoretical models. 

The SN provides us with two independent, yet equally 
useful means to reconstruct the mass-loss history. First, 
we are now observing a unique period in the forma- 
tion of SNR 1987 A as the high- velocity SN debris over- 
takes the slowly ex panding ER (Sonneb orn et alJ i 1998b: 
Sugerman et al. 2002). While this interaction will eventually 
destroy the circumstellar nebula, it provides the opportunity 
to explore the CSE by detailing the material directly through 
shock-heating. Furthermore, the increasing X-ray and UV 
flux from this interaction will reionize unseen portions of the 
nebul ar structure, revealing m any of them for the first time 
llLuo. McCrav. & Slavinll 994). 

Second, light echoes from the SN have illuminated a sig- 
nificant volume of the CSE, and can be used to constrain the 
three-dimensional (3-D) structure, density, and composition 
of the stellar outflows. Three general structures have been 
previously reported. 

Close circumstellar material: The inner few light years sur- 
rounding SN 1987 A have been studied between days 750- 
1469 (after the SN; 1987 Feb. 24.23) by CKH95. Echoes de- 
tected within 3" of the SN were mostly found to lie on an 
hourglass-shaped bipolar nebula, inclined 43° ± 5° south, and 
rotated 3° ±5° east of the line of sight, while a smaller set 
were also identified as roughly coplanar with the ER but at 
larger radii. 

Napol eon's Hat: Napoleon's H at (NH) was originally iden- 
tified by Wa mpler et alJ as an arc-like feature (also 
described as an archer's bow by Crotts & Kunkel 1991) loop- 
ing roughly 6" north of the SN. It has been variously inter - 
preted as a disk coplanar t o the ER llWang & Wampler 1992), 
a truncated double cone ( Podsiadlowski, Fabian, & Stevens 
or a parabolic bow shock at the terminal edge of 
a region evacuated by the m ain-sequence (MS) progenitor 
ftVang. Dvson. & Kahnll9"93h . 

Large radius contact discontinuity: Recognizing that red 
giants and supergiants typica l ly hav e slow, dense stellar 
winds, Chevalier & Emmering ( 1989) proposed that the in- 
terface between the RSG wind and surrounding medium (ho- 
mogenized and rarified by a prolong ed MS wind) would 
be a "snowplow" discontinuity ( Casto rTMcCrav. & WeaveJ 
119751) . As these are dense and cold, the interface should 
be both long-lived and have a significant dust component. 
Chevalier & Emmering ( 1989) predicted that a stellar-wind 
bubble in pressure equilibrium with the surrounding medium 
should occur at a radius of roughly 16 ly for reasonable RSG 
wind parameters. This correspon ds closely to echoes re - 
ported by Bond et al. (1989.1990). ICouch & Malir] i ll 9891) . 
and Crotts & Kunkel (1989,1991) roughly 9" from the SN, 
which are interpreted as the contact discontinuity (CD) be- 
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tween the RSG/MS winds. 

Until now, most modeling efforts have focused on the last 
~ 20, 000 years of evolution (roughly the kinematic timescale 
of the rings and evolutionary timescale of the BSG stage) 
while, in large part, avoiding the previous 0.3-1 Myr of the 
RSG phase, during which the NH and CD features formed. 
Observations of the full circumstellar environment provide a 
unique opportunity to detail the products of the progenitor's 
mass-loss. A more complete set of boundary conditions is 
necessary to properly investigate the SN system, and to dis- 
tinguish viable models from unrealistic scenarios. 

Using the technique of PSF-matched difference imaging 
on 16 years of optical imaging data, we report an observa- 
tional effort, at the highest sensitivity and resolution to date, 
to understand more fully the entire circumstellar environment 
(CSE) of SN 1987 A, and to recreate a more thorough history 
of the progenitor's mass-loss. This paper describes the de- 
tection and analysis of light echoes within 30" of the SN as 
observed in ground- based and HST i maging. This work is also 
summarized in Sugerman et al.N200 5). followed by a detailed 
analysis of how our results constrain formation models and 
the progenitor's evolution. 

The theory and application of light echoes are presented 
in |2j including a dust-scattering model which is used to 
constrain the density and composition of echoing material. 
The data are described in ^5] followed by the data-reduction 
pipeline (|4}. The light echoes, their three-dimensional posi- 
tions, and the tools necessary to analyze them, are presented 
in |3 These data are used to build geometric models of ma- 
terial within ~ 30 ly of the SN, the details of which are ex- 
plained in HJOD A complete picture of the CSE, including 
constraints on the density, chemical composition, and mass 
of all circumstellar material, is made in ^5] It is our hope 
that the results presented here will serve as a significantly- 
improved set of constraints for hydrodynamic models of the 
CSE formation, for stellar-evolution models of the progen- 
itor, and for addressing the general question of asymmetric 
and bipolar stellar outflows. 

2. LIGHT ECHOES: THEORY AND MODEL 

When a light pulse is scattered by dust into the line of sight, 
an observable echo is produced, provided the pulse is suffi- 
ciently luminous and the dust sufficiently dense. An echo 
observed a given time after the pulse must lie on the locus 
of points equidistant in total light travel from the source and 
observer, that is, an ellipsoid with known foci. This sim- 
ple geometry, shown in Figure |2] directly yields the three- 
dimensional (3-D) position of an echo, uncertain only by the 
assumed distance to the source. A complete discussion of 
light echoes, including single-scattering models and the ob- 
servability of echoes around a wide va riety of cataclysm ic and 
other variable stars, can be found in Sugerman (2003)). For 
completeness, we briefly summarize the salient points below. 

2.1. Light Echo Surface Brightness 

Consider the geometry shown in Figure|2j in which a source 
at distance D from the observer emits a pulse of light of dura- 
tion t. This light scatters off dust located at position r from 
the source, and reaches the observer a time t after the arrival of 
the unscattered light pulse. Since D is generally much greater 
than r, the echo depth z along the line-o f-sight can be approx- 
imated by the parabola ( Couderc 1939) 

„2 

(1) 




View from space 



FIG. 2. — Cartoon schematic of a scattered-light echo. Echoes lie on the 
locus of points equidistant in light travel from the source and observer-i.e., an 
ellipsoid-and appear as arcs or rings on the plane of the sky, with a one-to-one 
mapping between their 2-D image positions (right) and their 3-D locations in 
space (left). Note that the distances are not to scale: D is much larger than z, 
such that the segments from Earth to the scattering dust are nearly parallel to 
the line of sight. 



where p = rsind is the distance of the echo from the line- 
of-sight in the plane of the sky, and 9 is the scattering angle 
(provided D^> r). 

Treating the dust as a thin planar sheet of thickness Az, and 
adopting cylindrical coordinates, the surface brightness of an 
echo at a given wavelength A is given by 



B S c(X,t) = F(X)n H (r)- 



cAz 



(2) 



AirrpAp 

To facilitate computations, we adopt units of years and light 
years, making c = 1 . We define the integrated scattering func 
tion 



S(X,fj)= / Qsc(X,a)<j g $(p,,X,a)f(a)da 



(3) 



where a is the dust-grain radius, <j g = ira is the physical cross 
section of a grain, and 2sc is a grain's scattering efficiency. 
The number density of grains at position r with radii between 
(a , a + da) is given by 



nd(r,a)da = nn(r)f(a)da 



(4) 



where «h is the number density of H nuclei, and f(a) is the 
grain-size distribution function. The total flux from the light 
pul se incident on the dus t is F( X) = J F(X,t')dt'. We adopt 
the Henvev & Greenstein ( 1941) phase function $. Since the 
light pulse and dust have finite widths, the observed radial size 
Ap of the echo yields an estimate of its line-of-sight depth 



(5) 
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Both the observed flux F{X) and area element pAp diminish 
as D~ 2 , thus equation for surface brightness is distance 
independent, as expected. If the echo and source are close, 
they should suffer roughly the same extinction. This work is 
re stricted to single-s cattering events, and the reader is referred 
to IChevalieJ (U286) for a discussion of multiple-scatterings, 
which are only of significance when the optical depth of the 
scattering material > 0.3. 

2.2. Dust Properties 

Mathis, Rump L & Nordsieckl i ll 9771) found that a mixture 
of carbonaceous and silicate grains following a density dis- 
tribution given by f( a) oc q~ 3,5 adequately m odeled Galactic 
starlight ex tinction. | We ingartner & Draine (2001, hereafter 
WD01) and Li & Draine (2001) have calculated dust scatter- 
ing and absorption properties for graphite and "astronomical 
silicate," 7 from which WD01 have constructed a new density 

7 Tabulated valu es of dust p roperties are available at 

|http : //www. astro .princeton . edu/$\sim$draine/dust/dust .diel .htr 
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FIG. 3. — Dust grain density distributions for carbonaceous (solid lines) 
and silicate (dotted lines) dust from WD01. Thick lines are for the LMC 
model with Ry = 2.6, be = 1 X 10~ 5 , and thin lines are for the Galactic 
model with Ry = 3.1, be = 6 X 1CT 5 . For clarity, and to allow the reader 
to estimate the mass present in each size range, data are plotted in units of 
a 4 dn gl 

tion from Mathis et all 119771) . 
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FIG. 4. — Time-integrated spectrum of SN 1987A from 59 days before to 
42 days after maximum light. 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
Type 
F470 ; 


i i i i i i i i i i i i i i i i i i 
Type 
F612; 


, , i l , , , i l , , , , l , i , , 
Type S 
F688; 












:i i i i i i i i i i i i i i i i: 
Type M : 
F470- 


:i 1 1 i i 1 1 1 i 1 1 1 i 1 1 1: 
Type M: 
F612- 


:i 1 1 1 i 1 1 1 i 1 1 1 1 1 1 1: 
Type M: 
F688- 








l l l i l l l i l l l l l l l 
f\ Type LI 
: F470: 

p TT^rr^-r^ 


l l l l l l l i l l l i l l l 
t: Type LI 
: ''\ F612; 

: 1,,,,- 


l l i l l l i l l l l l i l 
T Type LI 
: F6B8; 







45 90 



45 90 135 



135 45 90 135 
Scattering Angle 8 

FIG. 5. — The average integrated scattering function S(X,fi), integrated 
over the total spectrum F(X) of SN 1987A (Fig.|4j and filter profiles (Fig. 
171 . Line style designates composition: solid lines show carbonaceous dust, 
and dotted lines show silicates. Columns are calculated for a given filter, and 
rows are calcuated for LMC dust with grain-size range types L, M, and S, as 
noted at upper right in each panel. This integral gives the energy scattered 
toward the observer in unit time by a unit volume of gas with unit density 
located unit distance from the source. 



nous near the SN, and that extinction is small within the echo- 
ing dust (implicit in the single-scattering approximation), flux 
from the SN and surrounding light echoes should be reddened 
by the same amount along the optical path between the LMC 
to the observer. This assumption allows the direct compari- 
son of the observed echo fluxes to the dust-scattering model, 
w itho ut the additional uncertainty of extinction corrections 
( SI4.6> . The spectrum in Figure |4]is therefore not extinction- 
corrected 9 . 



distribution f(a) to fit observed interstellar extinction. We re- 
fer the reader to their paper for the actual functional form. We 
adopt their LMC dust model for R v = 2.6, and f) C = lx 10" 5 , 
which is shown by thick curves in Figure [3] For compari- 
son, their preferred Galactic dust model (case A, Rv = 3.1, 
be = 6 x 10~ 5 ) is also shown. The r eddenin g ratio is de fined 
as R v =A V /E(B-V) dCardelli. Clavton & Mathili98"3) . and 
be is the total C abundance per H nucleus. 

2.3. The Light Pulse from SN 1987 A 

Usin g the early optical light curves from Hamuy et al. 
(1988), SN 1987Areached maximum light in V and R roughly 
87 days after its discovery on 1987 Feb 24. The light- 
curve duration r is measured as the effective width r = 
J F(X,t)dt /f(A,f m ax)- A homogenous time series of spectra 
from the first thousand days are publically available from the 
SUSPECT supernova archive 8 . We integrated these spectra 
over the first 420 days, and found an effective width of 100 
days in the wavelength ranges of interest. Integrating again 
over just this width, centered at maximum light, we gener- 
ated the time-integrated flux F(X) that will enter equation (0, 
shown in Figure [4] Note that at these early times, spectral 
coverage ended at 7300A, and there is no empirical spe ctru m 
available in the wavelength range of our reddest filter ( $3.21 . 

Referring again to Figure |2] since D 3> r, the trajectory 
from source to dust to earth is very similar to the direct source- 
earth path. Under the assumptions that the ISM is homoge- 

8 |http : //bruf ord. nhn . ou ■ edu/$\sim$suspect /index . html] 



2.4. The Light Echo Model 

To calculate the total echo surface brightness received in a 
given bandpass, equation must be integrated in wavelength 
over the filter response function T(X). The integrated scatter- 
ing function S(X,pi) must also be integrated over grain sizes, 
using equation Q and the dust-density models from WD01. 
Removing the geometric factors and gas density from equa- 
tion Q, we model echo surface brightness by computing 

1= JJ T(\)F(X)Q S c(\,a)a g <5>(n,\,a)f(a)dad\. (6) 

which is integrated numerically using Simpson's Rule 
JPress etalJl992l) . This gives the energy scattered toward the 
observer in unit time for unit volume of gas with unit density 
at unit distance from the source, and therefore it need only be 
scaled to the appropriate physical values to yield the predicted 
surface brightness. 

Figure[5]shows this integral for the three filters F470, F612 
and F688 (Fig.0. Results are given separately for carbona- 
ceous and silicate dust. Since the range of dust-grain sizes 
in the CSE are unknown, we integrate over three size in- 
tervals, denoted "S" (small, a = 3.5 x 10" 4 -0.01^m), "M" 
(medium, a = 3.5 x 10~ 4 -0.1/im), and "L" (large, a = 3.5 x 
10~ 4 -2.0/im). Th e type-L dust includes the full range of grain 
sizes in f(a) from Weing artner & Draind (120011) . 

9 Spectra of SN 1987A in the SUSPECT archive were dereddened assum- 
ing E(B — V) = 0.16 and Rv = 3.1, which we have removed to reproduce the 
observed values. 
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Scattering for grains with a < 0.1A/27T occurs within the 
Rayleigh limit, for which S oc A -4 . These small grains have a 
small scattering efficiency and scatter isotropically, while dis- 
tributions including large grains predominantly forward scat- 
ter. It should be noted that while S oc A" 4 for Rayleigh scat- 
tering, the largest color shift possible for the integrated SN 
spectrum (Fig.|4ji in our continuum filters is S oc A -3 , as can 
be easily calculated from the type-S dust in Figure[5] 

To apply this scattering model, colors of echoes are com- 
pared to the model values in Figure [5] to constrain the grain- 
size distribution, and the relative populations of carbonaceous 
and silicate grains. Dust (or gas) density is calculated from 
the scaling between the model and observed surface bright- 
ness, using the geometric factors p and Ap, measured from 
the data, with r=100 days, as determined above. 

Distinguishing between type S and L dust can be accom- 
plished a few ways. First, as noted above, small grains scat- 
ter less flux than larger ones, thereby requiring a higher gas 
density for a given observed surface brightness. Since greater 
density dust has stronger extinction, practical constraints exist 
for upper limits to gas density. If a single structure is probed 
at many scattering angles, and one assumes the density of that 
structure is constant, the variation of surface brightness with 
9 can also constrain the grain sizes. Dust composition can 
be constrained by comparing colors of observed and modeled 
echoes, although carbonaceous and silicate dust have nearly 
identical scattering properties at these optical wavelengths. A 
better probe of chemical composition are UV echoes, since 
carbon and silicate dust grains have strongly-differing fea- 
tures at A < 2500A ( Sugermanl2003l) . 

As a final note, in treating echoes with a planar approxi- 
mation, this scattering model assumes that the observed echo 
is well represented by a slab of dust of thickness Az, all of 
which is scattering the integrated flux F(X) toward the ob- 
server. In reality, the echo is the integral of many volume 
elements of gas, each of which scatters a parti cular time- 
depen dent flux from the SN lightcurve (see §2. 1 of ISugermanl 
2003). For echo observations in which a resolution element 
is much smaller than the outburs t duration cr (such as those 
from V838 Mon Bond et al. 2003), a much more detailed scat- 
tering model can deconvolve the spatially-resolved dust prop- 
erties using the time-dependent lightcurve. However, such 
treatment is unsuited to circumstellar echoes from SN 1987 A, 
since the image resolution is larger than 100 light days. 

3. DATA 
3.1. Observations 

SN 1987 A has been observed many times per year 
since day 375 (after core collapse) in a ground-based 
campaign to monitor the appearance and evolution of its 
light pulse as it illuminates circumstellar and interstel- 
lar material. Some dat a from this campaign have previ - 
ously been reported bv ICrotts. Kunkel & McCarthvl (119891) : 
[ Crofts & Heathcotel ill 9911) : ICrotts & Kunkell i 19911) : and 
Crotts et alJ dl995[) in studying the circums tellar environment, 
and in IXu. Crotts & Kunkel 119911 [1991 and IXu & Crottsl 
(1999) for interstellar echoes. 

For the present application, we require data with high 
signal-to-noise, good seeing (arcsecond or better), and suf- 
ficient resolution to resolve most crowded stars and detail in 
the narrow surface brightness features. Those epochs of imag- 
ing that match these criteria are listed in Tabled Data were 
primarily taken on the 2.5-m (100-inch) Dupont telescope at 
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FIG. 6. — Ground-based images of SN 1987A from day 1998. Each image 
is 1 30" on a side, in a standard orientation with north up, east left. The inset 
at bottom left is 10" on a side, (a) Narrow-band Ha image. The inset shows 
the extended nebula surrounding the SN, as well as stars 2 and 3. (b) F612 
continuum-only image of the same field as (a). 



the Las Campanas Observatory (LCO) at platescales of 0."163 
and 0."260 pix" 1 , but one epoch from the Cerro Tololo Inter- 
american Observatory (CTIO) 4-m telescope is also used. 

Since its launch in 1990, HST has been used at least yearly 
to observe the SN, largely through the efforts of the Supernova 
Intensive Survey (SInS) collaboration. We used archival Wide 
Field and Planetary Camera 2 (WFPC2) data taken between 
1994 and 2001 for the current application. Those epochs in 
which echoes were discovered are also listed in Table[T] 

3.2. Photometric Bands 

SN 1987 A is located in a region of the LMC with consid- 
erable nebular structure, as shown in an Ha image from day 
1998 in Figure |6ji. Furthermore, the three-ring nebula sur- 
rounding the SN is the strongest source of line emission in 
this region. These extended surface brightness features create 
significant sources of noise when searching for light echoes, 
especially in the inner few arcsec around the SN, since the 
flux from three-ring nebula changes with time. To minimize 
the contribution of nebular-line contamination, four specially- 
selected continuum bands were selected for the monitoring 
campaign, centered at 4700, 6120, 6880, and 8090A. The lab- 
oratory transmission curves are plotted in Figure These 
bandpasses were chosen to avoid the brighter emission lines 
of the diffuse LMC clouds (e.g., Ha, [N II], [O III], and [S II]), 
as well as strong emission lines from the SN ejecta. This has 
the added benefit of suppr essing flux from th e central SN at 
early times by 1-1.5 mags (Crotts et al .119951) . 

In keeping with the HST filter naming convention, we will 
refer to each continuum filter by its central wavelength in nm, 
preceded by an "F", such that the 6120A filter is F612. Note 
however that these differ from HST filter names since the lat- 
ter are appended with the width of the filter (i.e. F656N). The 
effectiveness of these filters at suppressing nebular-line emis- 
sion is shown in Figure [6j>, the F612 image from the same 
epoch as panel (a). There is no detectable signal from the 
large cirrus-like clouds or the inner circumstellar nebula. An 
echo results from the SN light pulse scattering off dust, and 
since its spectrum had a bright continuum over the duration of 
its outburst (Fig.|4j, these filters will not suppress echo signal. 

4. DATA REDUCTION 

Data were reduced in the following stages: (1) CCD pro- 
cessing, (2) homogenization, (3) pipeline pre-processing, and 
(4) PSF-matched difference imaging. Each stage is designed 
to be minimally-invasive, thereby maximally preserving faint 
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TABLE 1 
Light Echo Observations 



Epoch 


Day 


Telescope 


Platescale 


Filters 


1988 Dec 14 


659 


LCO-100 




0: 163 


B,V,R,I,Z 


1989 Mar 15 


750 


LCO-100 


Of 163 


6023,6067,6967 


1989 Dec 18 


1028 


LCO-100 


Of 163 


F470,F612,F688,F809 


1990 Feb 14 


1086 


LCO-100 


Of 163 


F470,F612,F688,F809 


1991 Mar 4 


1469 


CTIO-4 


Of 124 


F688 


1992 Jan 15 


1787 


LCO-100 


Of 260 


F612,F688,F809 


1992 Mar 20 


1851 


LCO-100 


Of 163 


F612 


1992 Aug 14 


1998 


LCO-100 


Of 260 


F612,F688 


1992 Nov 19 


2095 


LCO-100 


0"260 


F612,F809 


1 QQ^ Knv 99 

177J ZZ 


2462 


T cd i no 


\J. zou 


F470 Fft1 9 F^SS FSOQ 


1993 Dec 27 


2498 


LCO-100 


0"260 


F470,F612,F688,F809 


1994 Jan 6 


2508 


LCO-100 


0"260 


F470,F612,F688,F809 


1994 Apr 1 


2593 


LCO-100 


0"260 


F612 


1994 Sep 24 


2769 


//S77WFPC2 


0"045 


F555W,F675W,F814W 


1994 Nov 25 


2830 


LCO-100 


0"260 


F612,F688,F809 


1995 Jan 7 


2874 


LCO-100 


0"260 


F612,F688,F809 


1995 Mar 5 


2932 


//S77WFPC2 


0"045 


F555W,F675W,F814W 


1995 Mar 10 


2935 


LCO-100 


0"260 


F612,F688,F809 


1996 Jan 12 


3245 


LCO-100 


0"260 


F612,F688,F809 


1996 Feb 06 


3270 


//S77WFPC2 


()"045 


F555W,F675W,F814W 


1997 Jul 10 


3789 


//S77WFPC2 


()"045 


F555W,F675W,F814W 




4800 6000 



6200 6B00 7000 

Wavelength A (A) 



8000 8200 



FIG. 7. — Filter transmission curves for the narrow-band continuum filters 
used to observe SN 1987A. 




FIG. 8. — Differences of two photometrically scaled images, showing (a) 
the direct subtraction, and (b) the PSF-matched subtraction. 



structure and flux. We begin with a discussion of difference- 
imaging, followed by the four reduction stages listed above. 

4. 1 . Difference -Imaging Photometry 

Consider the problem of subtracting two photometrically- 
scaled images to study time variable phenomena. Figure |SJ? 
shows a typical direct arithmetic subtraction. Even under 
identical circumstances (same telescope, detector, airmass, 
etc.) the point-spread function (PSF) across the detector can 
change, caused by e.g. flexure of the mirrors, movement of the 
filter, or change in atmospheric seeing. As a result, residuals 
remain around every star in the field, leaving holes and haloes 





FIG. 9. — Two different PSFs of identical area are subtracted, yielding the 
residual in the right-hand panel. 



such as those seen here. 

This is easily understood by subtracting two Gaussians of 
differing peak and FWHM, but identical total area, as shown 
in Figurel9] The residual has a bright halo in the wings and 
a hole in the center, since an uneven amount of flux was sub- 
tracted from each point along the profile. This can be ad- 
dressed a few ways. Each image can be degraded, by smooth- 
ing with e.g., gaussian filters. This is an invasive procedure, 
smearing out structure and detail that might otherwise be stud- 
ied. Rather, an optimal technique is to remap the flux in one 
PSF to that of the other in a uniform and consistent way. 

Our technique is to seek a convolution kernel k such that the 
input star ; can be mapped onto the reference star r, or i = r*k. 
From the Convolution Theorem, the Fourier Transform (FT) 
of this mapping becomes 

FT(0 = FT(r) x FT(fc) (7) 

from which we isolate the convolution kernel 



£ = FT 



FT(/) 



FT(r) 



(8) 



iTomanev & Crottsl i ll 9961) implemented this method of 
PSF-matching in an external IRAF 10 package difimphot (DIF- 
ference IMaging PHOTometry), or iginally designed for use 
with the M31 microlensing project (Crotts 1992), and subse- 
quently updated for the MACHO microlensing project (e.g. 

10 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 
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lAlcock et all 19951) . In its simplest implementation, empirical 
or model PSFs are used with equation [8] to build the kernel, 
which is then convolved with the input image, and the refer- 
ence image is subtracted away. We have implemented further 
modifications to employ the robust PSF-fitting and building 
algorithms of daophot (Stetsonl ll987l) . as described in Ap- 
pendix |A] See Alard (2000) for an example of alternative 
methods of PSF-matching. 

The same data entering the arithmetic subtraction in Fig- 
ure |8j? are shown after PSF-matching in panel (b). Nearly 
all stars have been cleanly subtracted. Since Poisson noise a 
scales as the square-root of signal, the data are slightly nois- 
ier at the positions of subtracted stars than in the surrounding 
sky. Nonetheless, this represents an improvement in the back- 
ground noise by a factor of ~ 2 for faint stars, and > 5 for 
those of median brightness. The brightest stars with linear 
profiles are well-subtracted except in the cores, where Pois- 
son noise is highest. We can reliably detect extended surface 
brightness only 1 — la above the background sky, while for 
noiser regions (near bright stars or other variable sources) this 
threshold must be evaluated on an individual bases. 

PSF-matching can not completely subtract away stars with 
non-linear or saturated profiles for many reasons. (1) These 
stars do not have a uniform profile that can be described by a 
single function. (2) The PSF of such a star does not properly 
represent how flux is recorded on the chip. (3) Since these are 
the brightest objects in the field, even if the average profile 
is well-subtracted, the residuals will be largest. Thus these 
represent untreatable sources of error even in the best of sub- 
tractions. 

4.2. CCD Processing 

Ground-based data were reduced within IRAF in the stan- 
dard manner. Images were column biased, checked for resid- 
ual bias and, if appropriate, a full bias frame was subtraced. 
Images were then flat-fielded. Bad columns, rows or pixels 
were not cosmetically fixed; however, most images had bad 
data within 1-4 pixels of each edge, which were reset to zero 
value. Pipeline-calibrated WFPC2 images were used directly 
from the data archive. 

4.3. Homo denization 

For a given epoch, all data in a single filter must be com- 
bined into a single "stack" image. To perform difference- 
imaging, this data must also be geometrically-registered to 
a common reference frame. This latter step is critical since 
registration errors lead to large stellar subtraction residuals, 
a significant source of noise and error. We chose to regis- 
ter all ground-based data to the 0"163 pix" 1 platescale, for 
the following reasons. The primary interest of this work is 
to detect and resolve faint structure, thus data need to be of 
high-resolution. Faint structures become brighter when re- 
sampled onto a lower-resolution grid, but the ability to resolve 
these from the background (stars, subtraction residuals, cos- 
metic defects, statistical pixel noise) is diminished with the 
width of such features. Resampling high resolution data onto 
a lower-resolution grid also degrades the FWHM of stellar 
profiles, and PSF-matching works best with a well-sampled 
FWHM> 2 pixels. 

Since the images taken at 0."260 pix" 1 have high signal-to- 
noise and the stellar PSFs are well-sampled, the data can be 
interpolated onto a higher-resolution grid without great loss of 
information in fa int features. This has bee n rigorously tested 
and discussed bv lFruchter & Hookl (12001 within the context 



of "drizzling" data to higher resolution. Individual frames 
were all registered to a common orientation, then stacked. 
Since most frames are slightly offset from each other over 
a single epoch, this is effectively the same procedure as driz- 
zling observations with sub-pixel shifts. 

4.3.1. Geometric Registration 

One epoch at 0."163 pix" 1 was chosen, and the F612 im- 
age was geometrically -registered to orient the (x,y) axes with 
north and west, using an astrometric soluti on derived from 
stars in the UCAC-1 and 2 catalogs ( IZacharias etalJl2000t 
2004), which we determine to be accurate to within 0"03 of 
the Hipparcos astrometric reference frame in this region of the 
sky. This is the absolute reference orientation. All other fil- 
ters within that epoch were registered to this orientation, then 
all other epochs were registered to the corresponding filter in 
the reference orientation. For those epochs with non-standard 
filters, the nearest equivalent was used. 

Registration transformations for ground-based frames were 
computed for a third-order polynomial with half-cross terms, 
and resampling was accomplished using a bicubic interpo- 
lated spline. Since centroids are the dominant source of er- 
ror in these transformations, we employ the techniques of 
crowded-field PSF-fitting, for which exte nsive software ex- 
ists in the daophot package (Stetson 1987). An empirical PSF 
is built using actual stellar profiles from the frame, as out- 
lined in the daophot package. This PSF is then iteratively 
fit to groups of stars using a non-linear least-squares mini- 
mization, and centroids are output to a file. Avoiding stars 
with saturated or non-linear profiles, we find PSF-fitting has 
an average centroid uncertainty of less than 1 % of the FWHM. 
Matched pairs of stars in both input and reference images 
were made using starmatch, a highly effici ent impleme nta- 
tion of the triangle-matching algorithm of Groth (198a). In 
98% of the frames, geometric registration was accomplished 
with RMS errors less than 0. 1 pixel in both dimensions. 

4.3.2. Image Stacking 

All registered images in a given filter per epoch are com- 
bined into a single "stack." This begins with determining a 
list of reference stars, which are uncrowded, have a clean spa- 
tial profiles, high signal-to-noise, and constant flux. These 
lists are used to build empirical PSFs, and determine the pho- 
tometric scalings between images for combining and differ- 
encing. 

As a general algorithm, reference stars are identified from 
a group of images. Bright and unsaturated stars are selected 
by applying flux cut-offs to the stellar identification algorithm 
daofind, while crowding and clean profiles are addressed by 
the daophot software (e.g. pstsel), and custom software in 
the difimphot package. Aperture and PSF-fitted photome- 
try is performed on this star list for all images in the group. 
The photometry provides an initial guess of the photometric 
scaling between images, which is applied to the flux of each 
star per image. Stars that vary by more than a user-specified 
amount are rejected, and the procedure is repeated. Stars that 
vary within 3% of their median value are treated as constant, 
and constitute the reference-star list. 

PSF-fitted photometry is performed on each image in the 
stack using the reference star list. This photometry is used 
to determine an additive (sky) and multiplicative (scale) scal- 
ing to remove the median sky and photometrically normalize 
each image to the best-seeing, highest signal-to-noise image 



8 



Sugerman et al. 



in the stack. Data with unusually-poor seeing or large pho- 
tometric scalings can be rejected, and the remaining images 
are combined using a weighted average, where each image is 
weighted by 



scale 2 • sky ■ FWHM 2 ' 

renormalized such that Yl w i = !■ Bad pixels, from cosmic 
rays, warm or dead pixels, cosmetic defects, or large statisti- 
cal outliers, can be rejected using sigma-clipping if there are 
> 10 images in the stack, otherwise we use the median. We 
found through trial-and-error that when sigma-clipping rejec- 
tion could not be used, rejecting the highest and lowest pixel 
value for a given position in the data was sufficient to remove 
most bad pixels. 

4.3.3. The Homogenized Data 

The result of the above is a homogenous set of registered 
images, one per filter per epoch, which maximized the signal- 
to-noise while minimially -resampling the data. This last point 
is important, since each spatial resampling degrades the image 
quality, thereby smearing out faint echo signal, making them 
more difficult to detect. 

Finally, note in Figure |6] that Stars 2 and 3 are crowded 
with the SN and the immediately-surrounding nebula. Star 3 
is known to be variable, and will leave subtraction residuals 
in many difference images, depending on its flux at the time 
of observation. Since the SN is fading with time, it too will 
produce confusing residuals at early times, when it was still 
bright in the continuum. We therefore PSF-subtract the SN 
and Stars 2 & 3 using daophot. This was performed two ways. 
First, we allowed the subtraction algorithm (allstar) to locate 
its own centroids for the three sources. Second, we generated 
a geometric transformation between the ground-based frame 
and the WFPC2 F656N (Ha and [N II]) field, in which the 
centroids of these three sources can be measured to within 5 
mas, and fixed the centroids to these transformed values. 

WF PC2 data were registered and stacked as described in 
Suger man et alJ (120021) . Tiny Tim model PSFs were used to 
subtract Stars 2 and 3, as well as all other bright stars with de- 
tectable diffraction spikes. Note that in this and the following, 
ground-based and and WFPC2 data are never intercompared. 

4.4. Pipeline Preprocessing 

As explained in ^A] pre-processing requires the selection 
of reference stars, building of PSFs, and creation of refer- 
ence images. For both the ground-based and WFPC2 images, 
reference-star lists were constructed from those stars that are 
present and photometrically constant in each dataset. High 
signal-to-noi se an d clean PSFs were built using these lists as 
outlined in 3A.2I In the current application of difference- 
imaging, we are searching for extended surface-brightness 
features that vary in time. As such, we require a reference 
image ( 3A.1> that is free of all echo signal. Since echoes are 
present in each epoch in the ground-based imaging, an itera- 
tive procedure must be adopted to generate such a frame, as 
described below. 

All stacked images in a given filter were median-combined 
into a temporary reference stack, using zero-point offsets (for 
the sky) and photometric scalings derived from photometry 
of the reference-star list. Eac h stac ked image was then PSF- 
matched to and differenced ( 3A.3> by this temporary refer- 
ence stack. For a given epoch, the difference images from 



all filters were co-added into a single echo image. If a faint 
structure is truly an echo, it will appear as positive signal in 
all continuum filters, while background noise is a random dis- 
tribution about a fixed value. Artifacts, such as flat-fielding 
errors or ghosts, can either vary or remain constant between 
different filters, and therefore can present themselves as false 
signals in these images. 

Light echo candidates in each echo image were identi- 
fied using a custom-written implementation of the friends- 
of-friends algorithm. The median and standard deviation are 
computed within a moving box oimxm pixels, and any pixel 
whose value is more than Na above the median is flagged as a 
"friend." Subsequently, every pixel touching that friend that is 
more than na above the median is a friend of the friend. Every 
friend-of-a-friend is found in the image, with different groups 
merging when they touch. From a practical standpoint, this 
works best on block-medianed images, in which each pixel in 
a 4 x 4 moving box is replaced by the median value in that 
box. This smooths out a great deal of statistical shot noise, 
making echoes easier to detect above the background. In gen- 
eral, we found that N = 3 and n = 2 were sufficient. 

Friend-lists surrounding bright stars or more than 30" from 
the SN were rejected. The remaining groups were checked 
against those from the two adjacent epochs to ensure the sig- 
nals are real. It is highly unlikely that a bright echo will ap- 
pear in only one epoch, and not a few months prior or later, 
since such a structure would have to be extremely thin on the 
plane of the sky and along the line of sight, as well as being 
fortuitously placed to be illuminated at that particular epoch. 
Additionally, groups were checked to ensure they did not re- 
sult from ghosts from nearby bright stars; as there are many 
such stars near the SN, many such ghosts were found. Those 
lists that survived all cuts were turned into pixel masks for 
each image in that epoch. 

In some cases, consistent structures appeared as positive 
flux in early difference images and negative flux later, or vice 
versa, indicating that an echo was present in the reference im- 
age. In such cases, stacks were differenced from other in- 
dividual stacks well separated in time, and friend lists were 
generated from these difference images. 

All stacks in a given filter were then again median- 
combined, this time with the echo pixels in each image 
masked out. The above steps were repeated with these new 
reference images, and very few additional echoes remained. 
The third iteration generated reference stacks which did not 
produce any detectable echo residuals. We adopt these as our 
reference, echo-free images. The same technique was used 
for both the ground-based and WFPC2 images; however, in 
the latter case, echoes disappeared from the PC chip by 1997, 
after which considerable data exist. Thus generation of the 
WFPC2 echo-free images was much more straightforward. 

4.5. PSF-Matched Difference Imaging 

PSF-mat ching and difference imaging proceeded as ex- 
plained in S1A.3I For the ground-based data, difference im- 
ages were made between each epoch in each continuum filter 
and the corresponding reference image, for all three kinds of 
stacks: those in which no PSF stars were subtracted, those 
in which the SN and Stars 2-3 were subtracted with floating 
centroids, and those in which the centroids were fixed by their 
astrometry. 

In many cases, no matter how the PSF was created, it left 
circular residuals around star 2, and to a lesser degree around 
star 3. The cause is straightforward: the empirical PSF is 
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scaled to match the flux of star 2, as are its wings. Since the 
wings of any empirical PSF are noisy (due to the low signal- 
to-noise of the data that comprise them), this noise becomes 
significant when scaled to the flux of a bright or non-linear 
star, such as star 2. One may either replace the wings a 
posteriori by, e.g., an elliptical gaussian (as we do in PSF- 
matching) or attempt to minimize this effect by carefully con- 
structing the PSF. The former option is not built-in to daophot, 
and rather than write a revised version of this software, we 
opted for the latter. 

In subsequent discussions of structures found within 3" of 
the SN, a ground-based result is not deemed significant unless 
it appears in both PSF-subtracted difference images. Further- 
more, a result that lies within a noisy PSF-subtraction residual 
is not significant if we cannot demonstrate that the structure 
is real and not an artifact of the noisy wings. This is easy to 
test by examining the non-PSF-subtracted difference image, 
or by PSF-subtracting other bright stars and examining their 
residuals. 

For WFPC2 imaging, the Tiny Tim model PSFs have no 
noise in the wings, thus there is no need to consider alter- 
nate PSF-subtractions. All epochs of data were differenced 
from the corresponding reference image in the three filters of 
interest, F555W, F675W and F814W. 

4.6. Photometric Calibration 

To study the density and com position of the dust producing 
the observed light echoes ( 32. 1> . it is necessary to photometri- 
cally calibrate the continuum filters. Since fluxes determined 
from broad-band Johnson filters may be a poor fit to our nar- 
row continuum bands, we instead determine the calibration 
spectroscopically, as follows. 

Space Telescope Imaging Spectrograph (STIS) spectra in 
G430L and G750L (blue and red) were taken through the 
2"-wide slit in 2000 November as part of a Director's 
Discretionary HST all ocation to study hot spot evolution 
(Suge rman et alj [ 20021). This pointing in cluded stars 4. 7, 
34 and 35 from Walker & Suntzeff dl990h . respectively 5. "5, 
9."7, 17."0, and 16"0 from the SN. Each star is reasonably 
un-crowded on our ground-based frames, allowing for accu- 
rate PSF-fitted photometry with daophot. All four are close 
to the SN on the sky, and assuming they are similarly located 
within the LMC, all five sources should should suffer roughly 
the same extinction. 

Pipeline-calibrated spectral images were combined and 
one-dimensional spectra extracted using optimal techniques 
packaged within stsdas. Since stars were located at different 
positions orthogonal to the dispersion axis, their spectra suf- 
fer constant, positional-dependent offsets, which we corrected 
empirically by matching absorption features. 

The red and blue spectra were integrated over the filter re- 
sponse functions T(X), as shown in Figure to determine 
each star's observed fluence. The average of these fluences, 
compared to the observed counts (per second) of the stars in 
the reference images, give the calibrated flux of one count per 
second. The difference between the four stars per filter are 
consistent to within ~ 5%. 

Note that this calibrati on is for observed fluxe s, and has 
not been dereddened. Walker & Suntzeff ( 1990) find that 
the field surrounding SN 1987 A is consistent with a redden- 
ing of E(B-V) = 0.17, while IScuderi et al.1 (Tl996) determine 
E(B-V) = 0.19 from HST Faint Object Spectrograph spec- 
tra of star 2 . These two papers assum e Ry = 3.1 for the 
LMC, while Mi sselt. Clavton. & Gordonl (Tf999) give an av- 



erage value of Ry = 2.6 using ten sightlines. However, if 
one assumes the SN, surrounding dust, and these nearby stars 
al l suf fer the same extinction, then the dust-scattering model 
( 32. 1> can be applied without any reddening corrections. 

5. ECHO MEASUREMENT AND VISUALIZATION 

In this section, we present the PSF-matched difference im- 
ages, describe how light echoes were identified and measured, 
and discuss analysis and visualization tools that will be used 
in this work. 

5.1. Light Echo Data 

All epochs of imaging listed in Table ^ were pipeline- 
reduced, PSF-matched and differenced as explained in 50] A 
subset of the resulting difference images are shown in Figure 
fTOl Recall that the WFPC2 PC has 3.6 times higher reso- 
lution per pixel and therefore a smaller field-of-view. Since 
HST observations were rarely taken at the same orientation, 
the images are clipped with respect to ground-based data, and 
the clipped regions change with epoch. 

For all observations, images are shown in a standard north- 
up, east-left orientation, with major tick marks denoting an- 
gular distances of 10" each. The innermost echoes were par- 
ticularly bright up to day 1469, and are washed out at the 
brightness-stretch used to display the background and outer 
echoes. In such cases, a 12"-square inset centered on the SN 
is shown at the top right of each panel, scaled in brightness 
to better-resolve the innermost echo signal. All images are 
shown with a linear greyscale scaling, the limits of which dif- 
fer between images to better show the echoes in each particu- 
lar epoch. 

Panels (a) and (b) were taken through the V-band and a 
narrow continuum filter centered at 6067 A, respectively. In 
the first case, this image (day 659) was differenced from a 
sum of the F470 and F612 reference images, weighted to ap- 
proximately reproduce the throughput of the Johnson filter. 
This scheme worked well for removing most faint background 
sources. We compared this image to a difference made against 
another V-band observation from 1990 Nov. 5 (day 1171) 
and found the background was sufficiently flattened using the 
combination of narrow filters. The non-standard narrow con- 
tinuum image from day 750 was differenced against the F612 
image, the result of which is of equivalent quality to other 
F612 difference images. In both cases, it was preferable to 
use the echo-free reference images rather than the same filter 
at a later epoch, since this latter option introduces additional 
confusion from echoes present in the later data. 

The data in panel (d) were taken in a dithered pattern by 
imaging the SN on each of the four extreme corners of the de- 
tector. The resulting mosaic has different noise characteristics 
depending on the amount of overlapping data. 

A number of the images, e.g., panels (/), (g), and (i), show 
the limitations of PSF-sub tract ing Stars 2 and 3 prior to differ- 
encing, also described in 34.51 Star 2 is especially bright and 
often has significant flux in its wings at a radius larger than 
can be adequately modeled by an empirical PSF. This can be 
well seen in panel {g). Additionally, the wings of the model 
PSF are intrinsically noisy at large radii, and the bright flux 
from star 2 made even small-scale fluctuations large, produc- 
ing subtraction residuals within the model-PSF radius. This 
is most obvious in panel (/). 

There are man y noteworthy artifacts present in these data. 
As noted in 34.11 the Poisson noise is highest at the position of 
a bright star, and will therefore leave a subtraction residual at 
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FIG. 10. — Selected 50" X 50" difference images. North is up, east is left, and major ticks mark 10". The position of the SN, Stars 2 and 3 are marked with 
yellow stars. Echoes are marked by colored circles. The inset shows the central 12" at a different color stretch to resolve the innermost echoes. Colors are 
explained in SI5.2.2I (a) V image from day 659. (b) 6067 A image from day 750. (c) F612 image from day 1028. (d) F612 image from day 1086. (e) F612 image 
from day 1469. (J) F612 image from day 1787. 



its core. Such residuals can be seen throughout the data. Non- 
linear and saturated stars leave much larger residuals, both 
as noisy cores and diffuse wings, such as can be seen in the 
upper left corner or in the bright-star cluster ~ 15" WNW 
of the SN in ground-based images. Ghosts of bright stars, 
caused by internal reflections within the detector optics, are 
particularly noticeable in panels (g) and (/), but can be easily 



identified since they generally occur with the same positional 
offset from all bright stars. 

WFPC2 data suffer from a host of image artifacts that com- 
plicate echo detection. Foremost are the extended diffraction 
spikes from bright stars, which can extend over 6" (128 pix- 
els) from the source. These have been included in the Tiny Tim 
model PSF, but in practice can never be completely subtracted 
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FIG. 10.— continued, (g) F612 image from day 2095. (h) WFPC2 image from day 2769. (I) F612 image from day 2830. (j) F612 image from day 2874. (k) 
WFPC2 image from day 3270. (/) WFPC2 image from day 3789. 



away. Linear banding from charge transfer of bright stars is 
also omnipresent. Stray light from a bright star positioned 
within a well-defined avoidance region outside the PC field of 
view will cause broad arcs and thick bands of flux (fully dis- 
cussed in the WFPC2 ISR 95-06), which may be seen as dark 
regions along the left side of panel (h). A smaller, roughly- 
horizontal dark arc on the right side, roughly 6" south of the 
SN, is also a bright-star artifact. There may be much more 
echo signal in this image than we could reliably identify (next 



section), due to these considerable sources of noise. 

5.2. Light Echo Detection 

Echoes are present at a variety of positions and over a wide 
range of surface brightnesses. To measure its position, an 
echo may either be specified by the total locus of bright pix- 
els, or as the central position of the surface-brightness pro- 
file. However, quantifying the structure of extended, diffuse 
echoes, particularly in the wings of the brightest ones, is dif- 
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FIG. 11. — Sample radial profile from PA~ 50° on day 1028 (black points), 
fit by two and three moffats in panels (a) and (b), respectively. Each individ- 
ual function is shown with a different line style, and the sum of all functions 
is drawn in black. 



ficult, for a variety of reasons. 

Identifying all pixels associated with an echo is comprehen- 
sive in giving a better representation of the complete spatial 
region producing the echo. However, it is important that the 
echo positions can be visualized in a meaningful way. We 
experimented with describing echoes three-dimensionally as 
the complete locus of all p ixels (identified with our friends- 
of-friends algorithm, 34. 4> within the surface-brightness pro- 
file, and found that when many epochs were rendered in 3-D, 
the result was very difficult to analyze, since background data 
were completely obscured by foreground material. Since vi- 
sualization of the data is critical to interpreting the positional 
information contained in the echoes, we instead adopted a dis- 
cretized measurement by which the signal is described as a 
series of convolved profiles. 

Consider a bright, resolved echo with wings that extend to 
roughly three-times the FWHM of a PSF. This could result 
from a large, single distibution of dust, or from the convolu- 
tion of many discrete dust distributions, each with their own 
density. In the first case, if we identify the echo only at its 
center, we reject information contained in the wings, namely, 
that there is dust there. Rather, if we treat the extended flux as 
the convolution of many discrete echoes, even if there is only 
one echo present, we will better represent the spatial dust dis- 
tribution. 

The PSFs in most of our images are better described by 
moffats than Gaussians. The functional form for a moffat is 
given in Appendix lB.il it is essentially a modified Gaussian 
that can have a more rapid or shallow profile in the wings. 
Echoes were identified by fitting multiple moffats to radial 
surface-brightness profiles, generally measured by summing 
over wedges of arclength 10°. 

Moffats were fit using a highly "hands-on" implemen tation 
of Levenberg-Marquardt (LM) minimization (Appendix IB.U . 
Once the surface-brightness profile is measured, the user de- 
termines the radial range over which moffats will be fit, and 
can reject data that are known to be bad i.e., points coinci- 
dent with subtraction residuals. The user specifies the num- 
ber of moffats to be fit, and an initial guess for the peak, cen- 
ter, width, and exponent (3. Parameters are relaxed into the 
"best-fit" through iterative calls to the LM algorithm, with 
parameters replaced or reset when they deviate from reason- 
able values. "Reasonable" is fairly subjective but not without 
practical limits. For example, the width of a moffat was not 
allowed to be smaller than that of a stellar PSF, but was also 
not allowed to grow larger than ~ 2FWHM. 

Each profile is first fit with a single function, and additional 
moffats are added only as needed to reproduce the flux distri- 
bution. It is rare that more than three moffats were needed for 



a given echo. An example is shown in Figure^J The profile 
(points and error bars) is poorly described by a single Gaus- 
sian or moffat, as it has an extended plateau around p = 25 
pix, and will require at least two functions. The best-fit solu- 
tion for two moffats is shown in panel (a) in black, with the 
individual functions plotted in different line styles. Note that 
the fit is only marginal between p = 18-26 pixels. Subtract- 
ing these from the profile reveals a third peak centered near 
p = 18 pixels. Indeed, including a third function (panel b) is a 
much better fit, improving the \ 2 residual threefold. 

Once a successful fit is made, parameters are varied to en- 
sure the fit is robust. If so, the parameters are recorded, with 
the width of the echo Ap taken to be the moffat width a. Ra- 
dial profiles were examined at all position angles in each im- 
age, and out to a radius of at least 16" at early times, and 
32" at later times, to ensure that faint, distant echoes were 
not overlooked. The results have been marked in the differ- 
ence images, and color coded to compare echoes believed to 
originate in physical associations between epochs. 

As can be seen from the difference images in Figure fTol 
many echoes are faint, lie near bright stars, or in the case of the 
innermost echoes, near the positions of the SN a nd Stars 2- 
3. As with generating the echo-free images ( 34.4L faint echo 
signal was better resolved from the background by coadding 
the difference images from all filters in a given epoch. Echoes 
beyond ~ 5" of the SN were generally measured in these 
summed images. Identification of echoes within 5" is com- 
plicated by the unknown quality of the PSF-subtraction of 
the SN and Stars 2-3. Profiles were fit from difference im- 
age s ma de using both methods of PSF-subtraction explained 
in 34.31 and carefully compared for consistent results. When 
inconsistent, or we felt the subtraction of these sources could 
influence the result, the measurement was thrown out. In 
some cases echoes may be present but not reported. 

The difference image from day 659 (Fig. II Of a)) deserves 
special note. At this early date, the SN was sufficiently bright 
to saturate the detector long before sufficient surrounding sig- 
nal could be gathered. The PSF of the SN is too bright and 
extended to perform a PSF-subtraction, however we did find 
that subtracting a single moffat from the radial profiles left 
consistent echo signal in the PSF wings, as marked. 

5.2.1. Photometry 

Once identified as described above, echo brightness was 
measured as follows. All stars bright enough to leave a de- 
tectable subtraction residual were masked from each image. 
Echo fluxes were measured from these images with annular 
apertures of the same angular width as used during detection, 
and with radial widths given by the FWHM of each echo at 
each position. Since many echoes are extremely faint, we 
measured each twice, first as the surface brightness through 
the aperture allowing all data, and again counting only pos- 
itive pixel values. The corresponding median sky value was 
subtracted from each surface brightness, and the two values 
were rejected if less than the propagated noise, or averaged 
if greater. Counts were converted to fluence using the scaled 
exposure time from diffe renc e imaging and the photometric 
calibrations described in 34.61 

5.2.2. Discussion 

During their detection and measurement, echoes were cate- 
gorized as belonging to one of three expected structures, listed 
here in order of increasing radius on the plane of the sky: a 
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circum stellar hourglass-sh aped nebula reported in Crotts et all 
Jl995h . Napoleon's Hat (Wampl eret alJ 1 1990a). or circvT 
lar echoes at large radii (lBondetaLfTl989t ICouch & Malinl 
119891 ICrotts & Kunkelll 19891) interpreted as the contact dis- 
contin uity between the progenitor's RSG and MS winds 
(Chevalier & Emmering 1989). These have been color coded 
in Figure^] with the primary colors red, green, and blue, re- 
spectively. However in fitting moffat profiles, we found that 
additional identifications could be made between the hour- 
glass and Napoleon's Hat, generally colored cyan for echoes 
looping to the north of the SN, or yellow for positions south. 
Many of these southern echoes, once examined in 3-D, ap- 
peared at a larger distance, consistent with a southern coun- 
terpart to Napoleon's Hat, and have been marked in purple. 

That diffuse structure exists between the inner hourglass 
echoes (red points) and Napoleon 's hat (green p oints) has 
been previously reported. Crotts & Kunkel (1991) detected 
diffuse emission from 3" < p < 12" from the S N, measured 
in semi-circular apertures. ICrotts et"afl ill 9951) reported the 
detection of a "jet-like" feature, several arcsec in length, pro- 
jecting northeast from star 3 between days 1028 and 1727. 
These are consistent with the cyan and yellow echoes identi- 
fied in the same region of the current data, which now map out 
a much more extensive volume of the innermost circumstellar 
environment. 

Contact-discontinuity e choes h ave been p r evious ly re- 
ported by ICrotts et all (119891) . iBond et alJ {l990), and 
ICrotts & Kunkell l ) 19911) in observations only up to day 1028, 
while we are able to trace this structure through day 3270. 
Napoleo n's Hat was imaged as ear ly as days 850 and up to day 
1650 by Wang & Wampler ( 1992), while we detect its signal 
from day 659 to 2095. We believe this is the first discussion 
of the southern-counterpart (purple) to Napoleon's Hat. 

5.3. Visualization and Analysis Tools 

The 3-D analysis of these echoes is divided into the three 
broad categories defined above, and will be presented in |6j- 
[8] making use of visualization and analysis tools presented 
here. 

5.3.1. Coordinate Definitions 

In the following discussion of echo positions, we will make 
use of a few coordinate systems. On the plane of the sky, 
we use standard two-dimensional cartesian coordinates. The 
origin is always defined to be the SN, west is to the right and 
north is up. Polar coordinates p denote the radial distance 
from the origin, and position angle (P. A.) is measured as the 
counter-clockwise (CCW) angle from north. 

Three-dimensional coordinate systems are shown in Figure 
IT2l Let us define the standard orientation to be the orthonor- 
mal cartesian axes with z toward the observer, orthogonal to 
the 2-D image plane and centered at the SN, with x increasing 
west, and y increasing north. Thus, the line-of-sight distance, 
or depth, is measured by z. A point P in space can be de- 
scribed by its cartesian, spherical, or cylindrical coordinates. 
The spherical distance r is measured from the origin O, with 
polar angle <f> measured from north (y), and azimuthal angles 
are measured in the x— z plane CCW from the line-of-sight (z). 
We reserve an unscripted 9 for the scattering angle at point P, 

which is the angle between vector OP and z (under the reason- 
able assumption that d ^> r). cylindrical coordinates will be 
used for structures with rotational symmetry about an inclined 
axis z' . The point P has axial radius tt' measured orthogonally 




FIG. 12. — 3-D coordinate definitions for discussing light echoes. 



from z' ■ Transformed cartesian axes x' and y' may also be de- 
fined in the plane orthogonal to z' containing P. 

Rotational transformations are performed using a permuta- 
tion of the pitch-yaw-roll (PYR) convention, where rotations 
are performed about each of the three cartesian axes through 
the triplet (ip z ,9 x , (j> y ), where we first rotate (roll) about z by tp z , 
then about the rotated x (pitch) axis by 9 X , and finally about 
the rotated y axis by <j> y (yaw). Positive rotations are made 
about an axis following the right-hand rule. When referring 
to inclinations, the rotation is named for the axis about which 
the rotation is made, i.e., a rotation north or south is made 
about the x axis and named i x . The cylindrical axis z' is then 
inclined south/north by i x from the x - z plane, and rotated i y 
west/east from the y—z plane. 

5.3.2. Three-Dimensional Rendering 

To study echo positions in 3-D, we have custom-written 
graphics software to perform simple renderings, in which the 
positions of the echoes are transformed to allow viewing from 
any angle. An illustrative example is shown in Figure[T3l 

To realize a 2-D visualization of 3-D data, the following 
steps are performed. First, the user specifies the shift, rotation, 
scaling, and if desired, translated point of view of the observer 
and perspective. Using the matrix transformations explained 
in Appendix |Cj the data are transformed into a 2-D matrix 
of positions on the page, with the third dimension saved as a 
"z-buffer", which gives the depth of the point in or out of the 
page. Data are sorted by the z-buffer position so that points 
furthest from the observer are plotted first; in this way, distant 
points behind nearer ones will be covered. To facilitate an 
interpretation of depth, point size scales with distance from 
the viewer, with smallest points furthest away. 

In viewing 3-D structures, we find variation of point size is 
not sufficient to give the viewer an immediate sense of depth. 
We have added a simple ray-traced shading, assuming a sin- 
gle light source positioned behind the observer. In general, 
ray-tracing on an arbitrary shape is very difficult and involves 
complicated algorithms such as tesselation. However, all our 
echo data lie on parabolae of known functional form, thus 
the normal, or gradient, is always known. Points are shaded 
by the angle between the normal and line of sight, such that 
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FIG. 13. — Illustrative examples of 3-D renderings of light echoes in 
a simple bipolar nebula. This non-trivial example is intended to guide the 
reader in visualizing light echoes, and to annotate the rendering technique. 
The nebula, depicted at top, is inclined 45° to the line of sight, which is 
to the right; each lobe is 5 ly in radius, with the centers separated by 8 ly. 
Five echo parabolae are drawn, each occuring at the time (in years) indicated 
next to each curve. Below, renderings show the light echoes for the first two 
parabolae (first row), third parabola (second row), last tw o parab olae (third 
row), and all parabolae (last row). Axes are defined as in j!5.3.1l axis labels 
indicate the positive direction. The origin is set at the SN, and is always 
marked by a black dot. Major ticks mark 2 ly, and the origin is indicated by 
the longest tick along each axis. The coordinates and axes have been given 
a slight perspective transformation (Appendix |cj- Points are shaded using 
simple ray tracing, and larger points are closer to the observer. Left Column: 
Face-on view (the plane of the sky). Middle Column: Oblique view rotated 
45°. Right Column: Side view from far to the east. Note that the geometry 
of the nebula, as revealed in light echoes, is only clear when many echoes 
well-separated in time are considered at once. 




FIG. 14. — Regions of space surrounding SN 1987A probed by light echoes. 
The left and right panels show the view from far to t he no rth and east, respec- 
tively. The probable positions of the three rings ( 38.31 are centered at the 
origin. Circles of constant radii, in multiples of 5 ly, are drawn in light grey. 
The positions on the plane of the sky of Stars 2 and 3 are indicated by hor- 
izontal grey dashed lines. Echo parabolae corresponding to the epochs of 
observation used in this work (and listed in Table IT] are indicated in dark 
grey, with line style cycling from solid - dashed - dot-dashed - dotted - 
dot-dot-dot-dashed in order of increasing epoch. 



down. 

5.3.3. The Volume of Space Probed by Echoes 

Interpreting light echoes is often confusing, since they lie 
along paraboloids, a surface one rarely encounters when con- 
sidering geometric intersections. A plane, perpendicular to 
the line of sight, will appear as a concentric ring, growing in 
radius with time, until the paraboloid is larger than the mate- 
rial boundary. A thin, spherical shell of radius r, centered on 
the echo source, will appear first as a ring of increasing ra- 
dius as the paraboloids expand with time, until the paraboloid 
sweeps into the rear half of the sphere, after which the ring 
will collapse and eventually disappear at t = 2r. 

Other, more complicated geometries can appear in 
particularly-unrecognizable forms until the echoes have illu- 
minated a substantial part of the structure. For example, the 
inclined, bipolar nebula in Figure ^] can appear as (1) cen- 
tered or offset rings, both circular, elliptical, or highly dis- 
torted; (2) a distorted figure-eight; (3) multiple distorted rings; 
or (4) small curved segments, depending on the time and ori- 
entation at which echoes are observed. 

Figure El is more relevant to echoes from SN 1987 A, 
showing the positions of echo parabolae corresponding to all 
epochs of data in Table ^ Each panel shows a planar cross 
section of space passing through the origin: the x - z plane 
viewed from above on the left, and the y— z plane viewed from 
the side on the right. Concentric circles at radius intervals of 5 
ly are plotted to show how these echoes will illuminate struc- 
tures of different sizes. Also shown is th e probable geometry 
of the three rings, as determined in 38.31 and the line-of-sight 
positions of Stars 2 and 3 (as horizontal dashed lines), whose 
relative positions with respect to the supernova are unknown. 



0° (modulo 180°) is shaded light, and 90° (modulo 180°) is 
dark. In practical terms, that part of the paraboloid that is or- 
thogonal to the viewer's line of sight will appear whiter, and 
limbs will appear darker. For completeness, the line of sight 
in this context is orthogonal to the page, intersecting the ori- 
gin (black dot), and the observer is out of the page looking 



5.4. The Three-Dimensional Echo Positions 

Echo positions in 3-D are computed from the 2-D positions 
on the plane of the sky using equation Q. We have implicitly 
assumed that the echo center is produced by the light puls e 
maximum (which we take to be 87 days after discovery, §231 . 
and correct the observation date in Tabled by this amount 
when it enters equation Q. The angular size p is converted to 



The 3-D Map of SN 1987A 



15 



l|ii]|ili|iii|i.ii|lii|jiiIiiiJiil|iij|ili|iji|ui|iii|jii!iiiJiil|iij 



(a) 



««<Rs« **^!>& •»' 



I L I J IJ I.I LI Ll 1 I J I. 



(d) 







FIG. 15. — Rendered views of all light echoes identified in Fig. 1101 Points have been rendered using the method described in SI5.3.2I also see Fig. 1131 The 
left column shows observed views in the plane of the sky, the middle column shows views from the side (far to the east) and the right column shows views from 
the top (far to the north). Para bolae from the earliest and latest epochs at which echoes were observed are indicated in the side and top views. Point colors 
correspond to those in Fig. llOl except for the gold points in the top row, which denote echoes found in WFPC2 images which could not be resolved from the 
residuals of nearby bright stars in ground-based images. Note that the field-of-view changes between rows, but major tick marks always denote 2 ly. Top row: 
Contact discontinuity echoes. Middle row: Napoleon's Hat and southern counterpart. Bottom row: Inner circumstellar hourglass echoes. 



physical units assuming the SN of 50 kpc away. This value is 
a good average between many of the de rived distances to the 
LMC. which vary between 45-54 kpc (iGould & Uzalll998t 
lFeasil999tlRomaniello et al.l2000l) . 

The 3-D positions of all echoes identified in FigurelTUlare 
rendered in Figure^] using the methods described in §5.3.21 
The top, middle, and bottom rows show the data for the con- 
tact discontinuity (CD), Napoleon's Hat (NH), and the inner 
circumstellar (CS) material, respectively. Except for the gold- 
colored points in the top row (see below) colors correspond 
to those indicated in Figure ^3 A rigorous analysis of these 
data is presented in |j6H3 however we give a short qualitative 
discussion of the echoes below. 



5.4.1. The Contact Discontinuity 

In the plane of the sky (panel a), there is no data within p = 5 
ly of the SN. This results from the lack of echo imaging prior 
to 1988 December; this is seen in panels (b)-(c), for which 
there can be no data interior to the inner parabolae. There are 
two noticeable holes in the eastern blue echoes at (x,y) coor- 
dinates (-7,-1) and (-10,-6) ly, since a saturated star lies at 
each position. The bulk of the echoes are located between 20 
ly south and 4 ly north of the SN, with very few echoes posi- 
tioned further north. Since light echoes symmetrically probe 
all regions of space around the line of sight, this suggests that 
there is little structure in front and to the north of the SN, with 
the structure being illuminated inclined predominantly south 
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of the line of sight. 

The distribution of blue echoes east and west of the line 
of sight is fairly symmetric, however echoes extend to only 
p = 1 1 ly west, while they are present to the east up to p = 14 
ly. Comparing the two halves, there are eastern counterparts 
for all western echoes, yet there are no western counterparts 
bewteeen P.A. 250°-300° to the eastern echoes at p > 9 ly 
from P.A. 60°-110°. Since there is a cluster of bright and 
saturated stars in this western region, which makes reliable 
detection impossible, the absence of westerly echoes may not 
be meaningful. To the south, the echoes extend beyond p=14 
ly. These southern points do not appear symmetric about the 
y axis, but rather about an axis rotated 5°-10° eastward. 

When viewed from the side (panel b), the eastern/western 
points lie 4-14 ly in front of the SN, in a mainly vertical distri- 
bution extending 8 ly north and 12 ly south of the line of sight. 
The southern echoes are primarily oriented linearly along an 
inclined axis that points toward the SN. These echoes extend 
10 to 20 ly in front of the SN, corresponding to radial dis- 
tances of 14 to 28 ly. From this point of view, the southern 
echoes have a thin, sheet-like quality to them. 

Panel (c) shows the top view, from which we see that the 
eastern/western material follows a thin, curved surface that 
may be wrapping around the SN. These points lie at a fairly 
constant radius (in the x-z plane) of 12-14 ly from the SN. 
There is greater scatter in the eastern echoes than to the west, 
however both sets appear thin in the radial direction. This 
suggests the echoes may be the walls of a shell-like structure, 
such as a cylinder or spheroid. In contrast, the southern points 
again appear thin and sheet-like. 

The geometry of the complete structure is unclear from this 
subset of points. It may be a conical or prolate shell, although 
from only these data, it appears the structure may not wrap 
completely around the SN. Note that, aside from the small 
number of gold-colored points (see below) echoes do not ap- 
pear north of the SN. Instead, the material seems to stop 
abruptly along the northern edge of the early echo parabola 
(panel b). 

A small subset of the echoes have been colored gold. 
These points were identified in WFPC2 images very close 
to the bright northwest star cluster (P.A. 300°, p = 15"). As 
mentioned above, large residuals from these closely-spaced, 
bright stars made echo identification difficult-to-impossible in 
ground-based images. As they were observed at late times, 
these points lie 6-8 ly behind other echoes similarly located 
on the sky. No counterparts to these echoes are seen to the 
east, or at any other position angle at this radius and depth. 
These points do not immediately appear associated with the 
rest of the CD echoes, however we include them for com- 
pleteness. 

5.4.2. Napoleon's Hat 

Figures \15U-f show the 3-D positions of the NH echoes. 
As the northern echoes are already a known feature, we begin 
by describing them. These "NH-north echoes" are shown in 
green in Figure[T5l 

The NH-north echoes are reasonably symmetric between 
east and west on the plane of the sky (panel d). Viewed from 
the side (panel e), the bulk of the material lies along a thin 
shell that loops around the northern half of the three rings. At 
their closest to the observer, these echoes appear aligned with 
the equatorial plane containing the ER, and extend outwards 
roughly 1 .5 ly beyond the NOR. This shell looks to be aligned 
with the inclination of the rings, both south (as viewed from 



the east) and east (viewed from the north). From the side, the 
echoes could lie on a cylinder or spheroid aligned with the 
ring axis. Viewed face-on, the northern points are closer to 
the SN than to the east or west, suggesting that the structure 
has an elliptical cross section. 

Two nearly-horizontal features protrude from this shell to- 
ward the observeer, appearing horn-like when viewed from 
above (panel/). In the following discussion, we refer to these 
as the "horns," to distinguish them from the rest of the NH ma- 
terial. These are the echoes seen in early imaging that bridged 
between NH-north and the larger-radii CD echoes fFig.HOi. 

Between 1992-1996, arc-like echoes appear 1 .5 — 4 ly south 
of the SN. These are disguished from echoes arising in the 
same region at earlier times for two reasons. First, most of 
these echoes are isolated features on the sky, i.e., in radial 
profile they have distinct flux minima on either side of the 
peak, as opposed to the earlier southern echoes, which ap- 
peared as a convolved complex of nearby features. Second, 
the 3-D positions of these points are consistent with a thin 
shell-like counterpart to the NH-north echoes, while earlier 
echoes are widely dispersed around positions interior to this 
shell. 

We discuss these interior echoes within the context of cir- 
cumstellar material below, while the shell-like points are iden- 
tified as NH-south. These are shown in Figure ^] in pur- 
ple. Considered independently of the NH-north points, these 
echoes appear elliptical (panel d) and lie along a thin shell 
looping around the southern half of the three-rings. It there- 
fore appears that both northern and southern shells lie on the 
same structure, which (despite the temptation to call these 
"Napoleon's Collar") we will simply refer to as NH. 

5.4.3. Inner Circumstellar Material 

The 3-D positions of the CS echoes are rendered in Figures 
115b -/. Echoes passing within 1" of the SN wer e identified 
as part of the hourglass nebula (iCrotts et alJll995|) . and have 
been colored red. These form (nearly) complete loops at early 
times (up to day 1028), collapsing with time into increasingly- 
shorter arcs until only a single point was distinctly-identifiable 
on day 2462. Outside of these innermost loops, those echoes 
that form (nearly) complete loops or pass north of the SN have 
been color coded cyan, and those that remain primarily south 
of the SN are shaded gold. This particular choice of cate- 
gorization is arbitrary, and has been made only for clarity in 
viewing the echo positions. 

Viewed from the side (panel h), the red echoes do indeed 
appear to trace out a hollow, conical or tube-like surface. Vi- 
sually, these points appear inclined roughly 45° to the line 
of sight, with the opening oriented north and away from the 
observer, while viewed from the top (panel z), they appear to 
lie along an axis rotated 10-15° from z. That there are no 
loop-like counterparts in front and slightly north of the SN 
is a selection effect of the positions of the echo paraboloids, 
as shown in panel (h) and Figure ^] This latter figure does 
suggest that the earliest echoes may have skirted the south- 
ernmost edge of the hourglass. Panel (h) does show cyan 
and gold echoes positioned roughly at these expected posi- 
tions. There is, however, no a priori means to establish which 
of such echoes intersect the hourglass, versus which echoes 
are from material external to that structure. The gold echoes, 
when viewed from the side, are suggestive of a larger double- 
lobed structure surrounding the inner hourglass. Again, phys- 
ical association of these points is difficult to establish, and we 
defer further discussions to our quantitative analyses in ^8] 
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(a) (b) 

FIG. 16. — As Fig. ll5fc -/>, showing only the CD echoes, where now "shell" 
echoes are dark grey, "radial" echoes medium grey, and lig ht-grey points 
denote the WFPC2 echoes th at are difficult to interpret ( i|6.1.H . The probable 
positions of the three rings f i|8.3l are indicated in each panel. 
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FIG. 18. — CD light echo positions in spherical coordinates <j> versus radius 
r. Shell (dark grey) echoes are plotted with circles, filled for those east of the 
line of sight, and open for those to the west. Radial (medium grey) echoes are 
plotted as X , and the small subset of WFPC2 echoes to the northwest (light 
grey) by asterisks. The best-fit lines through the east and west shell echoes 
are indicated by solid and dotted lines, respectively. 



(a) (b) 

FIG. 17. — Schematic intersection of (a) outward-facing and (b) inward- 
facing cones with light echoes. The source is marked with a solid dot. 



6. CONTACT DISCONTINUITY ECHOES 

As explained briefly in ^ a strong body of evidence sug- 
gests that the progenitor star passed through a RSG phase be- 
fore terminating as a BSG. Those echoes identified as part 
of the contact discontinuity (CD) between the RSG and MS 
winds were color coded blue in Figure [K)] and blue or gold 
in Figure [T3b-c. Figure [Tol shows the same 3-D echo posi- 
tions, but in greyscale. In §5.4.11 we noted that the geometry 
of the eastern/western points appeared different from those 
to the south (P.A. 150°-210°). Here, the former points have 
been shaded are darker tone than the latter, and the WFPC2 
data for which there are no ground-based analogs are shaded 
the lightest. 

6.1. Geometry of the CD 

The reader may wish to skip to §6.2l for a summary. From 
the side (Fig. I16t>). the CD echoes appear to lie along a 
conical surface. A cone, opening away from the SN, and 
inclined roughly 10° south, will contain many of the echo 
points. However, the fit is only good for the southern mate- 
rial (medium-grey points). Consider the intersections of echo 
parabolae with cones, as shown in Figure^] Recall that the 
parabolae become larger with time. Light echoes through an 
outward-facing cone (panel a) intersect it at increasingly-large 
distance from the source with time, consistent with the south- 
ern points (medium grey points). However, eastern/western 
material (dark-grey points) remain at constant or decreasing 
distance with time, more consistent with a spheroidal shell, 
or an inward-oriented cone (panel b). This implied "saddle" 
geometry is unlikely, since no counterparts to the southern 
echoes are observed to the north. Qualitatively, the structure 
of which these echoes are a subset is unclear, and we proceed 
to a more detailed analysis. 

6.1.1. Geometry in Spherical Coordinates 



As the eastern/western echoes appear to form a curved sur- 
face surrounding the SN, let us first look at these data in the 
spherical coordinates <fi and r, shown in Figure^] Medium- 
grey points in Figures ^] are plotted as crosses, light-grey 
with asterisks, and dark-grey as circles, filled or open if they 
lie east or west of the line of sight. Roughly speaking, the 
circles lie in a horizontal band between 10 and 18 ly from the 
SN, while crosses are distributed vertically at constant <f>. A 
radial line will appear vertical in such a plot, while a sphere 
will appear horizontal. In the broadest terms, the circles seem 
to form a shell-like feature at roughly constant radius from the 
SN, while the crosses form a radial feature at roughly constant 
inclination. In the following discussions, we will distinguish 
these sets of points by referring to the crosses as a "radial" 
feature and the circles as a "shell." 

The radial echoes extend 14-28 ly from the SN. Their dis- 
tribution is nearly vertical about <j> = 130°, consistent with a 
radial line inclined 40° south of the line of sight. A few points 
are located at r < 17 ly and <f> = 120°. Whether these points 
are part of the shell, or yet another structure inclined 30° to 
the line of sight, will be addressed in the next subsection. 

To the west (open circles), the shell is radially thin at all an- 
gles. At any given <$>, these points show a scatter of roughly 2 
ly, which combined with the intrinsic echo width (Az ~ 1 — 3 
ly; see Sugerman 2003) suggests a shell thickness of ~ 4 ly. 
To the east (filled circles), the shell thickness increases from 
a scatter of 1 ly at large </> to up to 6 ly at small <f>. Both 
eastern and western sets of points are visually suggestive of 
curved structures, nearest to the SN around <f> = 90° (i.e., the 
x—z plane), and flaring outward in radius above and below 
that plane. To test this, we first examine how well each set is 
described by a linear fit. These are shown as solid and dotted 
lines for the eastern and western points, with an RMS scat- 
ter about each line of 1.1 and 0.8 ly, respectively. Both lines 
are generally a good fit to the data, however at <fi < 90° both 
sets of points deviate to larger r. Points at <fr > 120° may 
marginally deviate from these lines to larger radii. 

The western half of the shell is closer to the SN by 2.7 ± 1 .4 
ly, using the best-fit lines described above. Although only a 
2-a difference, this could indicate an asymmetry in this struc- 
ture. If the CD, whatever its actual shape, has symmetry, then 
its axis is spatially offset from the SN. If we require this axis 
to pass through the SN, then the asymmetry may lie in the hy- 
drodynamic interaction that formed this structure, including 
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the outflows from the progenitor, and the medium into which 
they propagated (e.g. interaction with the remnant of a previ- 
ous, nearby SN). 

A cartesian line that does not intersect the origin will ap- 
pear as a convex curve in spherical coordinates, symmetric in 
slope about the minimum in r. That the shell points appear 
to lie on such curves suggests they are part of a nearly-linear 
surface of revolution, such as a cylinder or a cone. Comparing 
Figure [H5]to the expected distributions for various geometric 
surfaces, shell points could be consistent with an inclined el- 
lipsoid, cylinder, or inward-opening cone. Heuristically, the 
inward-facing cone may be the best fit since its apex translates 
into a highly inclined, nearly vertical distribution, as seen for 
the radial points. 

The northwestern WFPC2 echoes, marked as stars in Figure 
[TSl are nearly colinear with the western shell material. This 
could be consistent with the conical distribution described by 
the rest of the CD echoes if the cone becomes more cylindri- 
cal as it passes the SN. However, it is difficult to reconcile 
this structure with the shell material that flares to larger radius 
at cf> < 90°, indicative of a cone growing to larger radii as it 
nears the SN. The association of these points with the rest of 
the CD echoes is unclear, and they may instead arise from iso- 
lated inhomogeneities in the surrounding ISM, or an unusual, 
extended image artifact, rather than from stellar outflows. 

6.1.2. Finding the Orientation 

If the shell and radial material are part of a spheroidal or 
conical structure, it should have a well-defined axis of sym- 
metry. We assume this axis passes through the SN, and at- 
tempt to find it by \ 2 minimization. Note that only a portion 
of the physical CD structure has been illuminated. A numeric 
minimization of some merit function will always be more suc- 
cessful at recovering the axis of an object with a complete 
cross section (e.g., a full cylinder), than with a partial one 
(e.g., a cylindrical segment). Thus, exploring this material's 
geometry is inherently uncertain. 

First, we directly fit the axis by finding the line that min- 
imizes the geometric distance to all echo points (Appendix 
IB.3i . Constraining the fit to have zero offsets, we find the 
best-fit line is inclined 37° south and rotated 8° east of the 
line of sight. The formal error of this fit is 0.°5 per axis of ro- 
tation, and while it is robust to perterbations in the parameters, 
we expect the actual uncertainty is higher. 

The echo points are transformed onto a cylindrical coordi- 
nate system about this line, designated by the primed coor- 
dinates in Figure [21 Figure^] shows the "radial profile" of 
axi al po sition versus radius for the transformed data. As noted 
in 35. 21 we only measure the center of echoes, therefore points 
have 1—4 ly incertainty in their spatial extent. Symbols are the 
same as in Figure^] except for the "+" marks, which we dis- 
cuss shortly. In these coordinates, a cylinder will appear as 
a horizontal feature, a regular cone as an inclined line, and a 
biaxial ellipsoid (with major axis along z') as a curved feature 
with dir'/dz' < and d 1 n'/dz' 1 < 0. We have again plotted 
the best-fit lines through the shell material, which show that 
the bulk of it lies on a linear incline, and thus delineates a 
conical structure. 

The shell material is best fit by the line n' = (14.4 ±0.4)- 
(0.5 ± 0.04)z', giving the conical structure an opening angle 
of ~ 53°. Considered independently, the eastern and west- 
ern halves are fit with zeropoints of 16.2 ± 0.5 and 13.3 ± 0.5 
ly, and slopes of -0.62 ±0.05 and -0.43 ±0.05, respectively. 
Although these values differ, the eastern points appear only 
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FIG. 19. — CD echoes, transformed into cylindrical coordinates with z' 
along the axis of symmetry given by i x = 37° and i y = 8°. Coordinates are 
shown in Fig. 1121 with z' along the axis, an d it 1 measured radially from it. 
Symbols and line-styles are the same as Fig. 1 181 with the heavy solid line is 
the best fit through all data. A subset of the radial echoes have been marked 
by +, indicating "cylindrical jet" material. 



marginally offset to higher radii than those to the west. Close 
to the SN, the shell flares outward to a larger radius of 15-16 
ly, as first suggested in Figure ^] The cone continues with 
roughly constant slope and thickness until shell data ends at 
z' = 16 ly. Much of the radial material is colinear with the 
shell echoes between z' = 15 and 28 ly, but with slope -0.25, 
for an opening angle of ~ 30°. Treating the shell and radial 
points as part of a single "conical-shell" structure, we find that 
it does have a curvature to its profile. 

A small subset of radial material between z' = 1 1 —26 ly lies 
at 7r' just below the sh ell ma terial. After considering many 
possible inclinations ( SI6.1.3> . we have identified the points 
marked by "+" symbols as behaving differently from the bulk 
of the radial data. If we consider these points to be part of 
the conical shell, then the entire ensemble of radial points ( x 
and +) have the same slope as the shell-only echoes (circles). 
In such a case, the conical shell thickens with distance from 
the SN, and approaches an apex or termination at a distance 
of r = 28 ly. The subset of radial points in Figure^Dthat are 
consistent with a 30° inclination ( 36.1.U are all found within 
this new set of points. If we treat them as an independent 
feature, they are well fit by a horizontal line. These echoes 
lie within tt' < 3 ly, and extend inward toward the SN, sug- 
gestive of a quasi-cylindrical region. If a narrow jet pierced 
through the CD, its interaction with the RSG wind might have 
formed such a narrow conical or cylindrical high-density fea- 
ture. With this in mind, we will treat these points as a separate 
feature from the conical shell, and will refer to them as the 
"jet- 
Finally, the WFPC2 echoes northwest of the SN are located 
at a radius inconsistent with the slope of the other CD echoes. 
Since there are no echoes located between z' = and 4 ly, 
one can only speculate how the profile of the CD structure 
changes as it passes the equatorial plane. If, as noted above, 
the shell flares out to large radius at small z' , these WFPC2 
echoes do not lie on the boundary of a structure, but instead 
may simply be part of a local overdensity of gas and dust, as 
noted previously. If the profile does curve in, and it contains 
these echoes, then one must ask why no material joining the 
shell to these points was illuminated. These echoes appear 
isolated, and have no counterparts at other P. A.s, thus we give 
them little weight when considering the CD structure. 

Since the CD echoes appear to lie along a conical shape, 
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its inclination should be better studied by performing the \ 2 
minimization of a cone rather than a single line. Unfortu- 
nately, the minimization of this merit function (Appendix lB.4> 
provided few useful results. Because the echoes trace out 
only a limited portion of the structure on which they are lo- 
cated, the best-fit cones have inclinations larger than i x = 60°. 
Such highly-inclined cones are poor fits to the total dataset, 
but minimize \ 2 by passing very close to a large subset of 
points at large distance from the SN. As a result, the remain- 
ing discrepant points effect the merit function very little, even 
though the overall fit is poor. While these fits give little infor- 
mation on the inclination of the conical structure, inclinations 
in the range of i x = 30 - 50° consistently require rotations of 
i y = 5-10°. In the following subsections, we describe alter- 
nate techniques for establishing the most probable orientation 
of this conical structure. 

6.1.3. Alternate Orientations 

As it cannot be established whether the proposed inclina- 
tion is correct, we adopt a different approach by studying how 
the radial profiles vary over a range of inclination angles. We 
consider inclinations i x between 30° and 56° from the line of 
sight, tested every 2°, with the corresponding rotation i y given 
by the best-fit 3-D line to the data at each inclination. For 
clarity, the CD echoes are rendered for i x = 30°, 40°, and 50° 
in Figure [20] The jet material has been shaded light grey, to 
distinguish it from the radial echoes. Qualitatively, the echoes 
appear to lie on a conical surface of revolution only for the 40° 
and 50° inclinations. An inclination near 40° has the most 
overlap (in projection) of the jet echoes and the SN, which 
strengthens its intepretation as a narrow, cylindrical jet. At 
50°, the radial and jet material lie neatly on a surface of cir- 
cular cross-section, suggestive of demarcating the edge of an 
evacuated, cylindrical region, perhaps connecting to the ER 
and SOR. Although this inclination is not suggested by the 
previous analysis, we question whether it could be a viable fit 
to the full set of CD echoes. 

As in Figure^] we transform the CD material into cylindri- 
cal coordinates about each inclination axis, and generate ra- 
dial profiles through which we measure the best fit lines. The 
slopes of those lines are summarized in Figure |^ The slope 
(dir' /dz') of the shell material is fairly uniform between -0.45 
and -0.51, with this steepest slope occuring at i x =42° -44°. 
Although the variations appear smaller than the quoted error 
bars, we note that uniform positional errors of unity were as- 
sumed in the least-squares fits. The radial echoes are nearly 
horizontal at the smallest and largest inclinations, with the 
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FIG. 21. — Slopes of the best-fit lines 7r' = a + bz' through the CD echoes, 
inclined to the south as given along the abscissa. Symbols denote the same 
material as Fig. 1191 except all shell material is represented by open circles. 
Grey error bars show the formal uncertainty from the least-squares fit, assum- 
ing equal errors for all data. 



steepest slope also occuring near i x = 42°. At no point are the 
slopes of the radial and shell e choes equal. This reinforces 
the interpretation from 36.1.21 that radial points, if axisym- 
metric, lie on a structure with much shallower opening angle 
than that of the shell material. The slope of the jet material 
is very sensitive to inclination, with the implied axisymmetric 
structure appearing as a wide cone opening away from the SN 
at i x = 30°, a narrow cylinder at i x = 36°, or a narrow cone 
opening toward the SN at i x > 38°. Are these results consis- 
tent with viewing the hypothesized conical shell at different 
inclinations? 

Consider viewing a complete frustum at different inclina- 
tions, as sketched in Figure|22] The radial profile of tt' versus 
z' is linear only close to the true axis, marked as line A. Away 
from A (e.g., line B) the plots become highly non-linear, and 
the slope of the best-fit line becomes more vertical, with in- 
creasing inclination away from the true axis. This is shown in 
Figure 123b . which is generated by viewing a simulated frus- 
tum at different inclinations. Here, the slopes are steepest at 
the highest inclinations, contrary to the data in Figure|^ Ex- 
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FIG. 22. — Cartoon sketch of a frustum, and its intersection (/,/') with the 
horizontal plane H. The axis of symmetry is line A, common to planes V 
and H, where V and H are orthogonal. We consider inclinations i of the axis 
relative to line B, also in plane V. 
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FIG. 23. — Data similar to Fig. l21l but for a simulated frustum, oriented as 
in Fig. 1221 with its axis along line A, and inclinations i x made in the plane 
V. Panels (a)-(c) use a frustum with large opening angle, simulating the shell 
echoes. Best-fit slopes are calculated as follows: (a) Points sampling the 
complete surface of revolution, (b) Points within 30° of plane H. (c) Points 
from 30° below to 60° above plane H, meant to simulate the positions of 
actual shell echoes, (d) Points along the upper-half of a frustum with small 
opening angle, simulating the jet material. 



amination of the renderings in Figure [2U] shows that for the 
shell and radial echoes, material is primarily located left and 
right of the symmetry axis, not above and below. Since the 
inclinations move this axis through the region with no echoes, 
there is no material present to cause the effect just described. 

Rather, the situation is akin to rotations in the plane of re- 
flection between two curved sheets. Reconsider Figure l22l 
in which two lines I and /', with refiectional symmetry about 
plane V, lie in plane H. The slope of line / as measured from 
line B (inclined in V) will be similar to, but shallower than, 
the slope measured from line A, for moderate inclinations. If 
we instead consider thin, curved sheets located on the frus- 
tum (similar to the shell and radial material illuminated by 
echoes), but only extending some angle above and below H, 



the same effect will occur. This is shown in Figure l23t>. Un- 
der these circumstances, the slopes are the shallowest at the 
highest inclinations, consistent with the conical-shell data. 

If the shell and radial material were unformly distributed 
above and below the axis of symmetry, then Figure I2T1 would 
imply the actual inclination is near 43°. However, as seen 
in Figure |20| more material is above the axis than below it. 
Figure|23j: shows the expected variation for a frustrum wedge 
chosen to approximate the actual distribution of data. Now 
the steepest slope occurs at an inclination a few degrees larger 
than that of the true axis, leading us to favor an inclination 
~2°-A° less than 43°, i.e., near i x ~ 40°, for the conical-shell 
material. 

If we interpret the jet material as a segment of narrow, semi- 
cylindrical sheet, its change in slope should appear as Figure 
1231 /. Similar in appearance to the jet data in Figure|^ these 
simulated points also suggest the extremal slope is offset by 
2°-4° from the true axis. If we continue to interpret the jet 
as an independent feature, its inclination may also be in the 
range 38°^10 o . 

6.1.4. The Conical-Shell Cross Section 

From the previous arguments, the favored inclination of the 
conical shell is around i x = 40°, although the alternate val- 
ues of 37° and 50° have yet to be ruled out. If the echo 
points do lie on a surface of revolution, the cross section as 
viewed along its axis should appear circular. An elliptical 
cross section would either indicating a viewpoint inclined to 
the symmetry axis, or that the surface is inherently elliptical. 
We investigate this by fitting ellipses to cross-sectional pro- 
files through x 2 minimization, as explained in Appendix lB.21 
Since the axial radius ir' varies with z' for the conical-shell, 
points must be normalized to a common radius. The radial 
profile at a given inclination is divided into bins of width 1 ly 
in z, then axial radii are normalized by dividing out the aver- 
age value in that bin. 

For each inclination tested in the previous subsection, the 
cross-section of the e choes is calculated (viewed along the 
axis) as described in 36.1.21 A subset of these are shown in 
Figure [24] for the inclinations noted in each panel. The shell 
echoes lie closely scattered along two arcs of a circle or el- 
lipse. With increasing inclination, less material appears under 
the abscissa. At small i x , the radial points appear to form the 
bottom edge of the ellipse, and become intermixed with the 
shell echoes at higher inclination. The jet material form a 
smaller circular shape at moderate inclination, and eventually 
merge with the other material into a single ellipse at i x = 56°, 
supporting the suggestion from Figure l2T)l c) that the inclina- 
tion axis is large. 

Assuming one of these to be the true inclination, we can 
build the conical shell implied by its radial profile, then test 
whether the other cross sections are consistent with inclined 
views of that data. This is shown in Figure [25] Linear fits 
to the radial profiles at inclinations of 37°, 42°, and 50° de- 
fine the function tt'(z') for the shell and radial data. For each 
inclination, the function is rotated about the z' axis through 
position angles chosen to approximate the actual cross sec- 
tion. Points along these surfaces are chosen at random, and 
given random scatter, to form a simulated conical shell. Each 
shell is shown in a single row, as labeled on the left side, and 
inclined north or south of z' by the angle given at the top of 
each column. 

Suppose the CD is actually inclined at i x = 37°. That panel 
in Figure |24] should be compared with the 0° inclination of 
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-1 o 1 
x' (normalized units) 

FIG. 24. — Cross sections of CD material, viewed along a hypothetical symmetry axis with inclination as no ted in each panel. All data are norm alized by their 
average distance (ir 1 ) measured in bins along z' ■ Labels x' and y 1 are the cartesian axes denned in Fig. 1121 Symbols are the same as Fig. |WI except all shell 
material is designated by open circles, which are s haded grey to facilitate the viewing of other symbols. Comparison of this figure with Fig. l25l shows the most 
probable CD inclination is near ;' v = 40° (see text, i|oT4l 
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FIG. 25. — Cross-sectional views similar to Figure l24l but for simulated conical shells. Each shell reproduces the CD data viewed at the incinations listed at the 
left of each row, and is viewed at the inclinations listed at the top of each column. Open circles represent simulated shell material, and simulated radial material 
is marked by X. Also plotted are the best-fit, 5-parameter ellipses and centers through the full dataset (thick lines, solid star) and only through the shell points 
(thin lines, open star). The ellipticity bja for each fit is given at top right and left of each panel, respectively. 



the top row in Figure [25] it is a reasonable fit since it was so 
constructed. When inclined north or south, the data should 
transform as the simulation does. If inclined to the north, the 
data in Figure [24] at i x = 30° is reasonably similar to the sim- 
ulation at 6.°5 in the top row of Figure[23] When inclined to 
the south, however, the simulated radial material in Figurel25l 
appears above the abscissa with considerable scatter, a poor 
fit to the data in Figure [24] at i x > 43°. Thus, an inclination 
of 37° seems unlikely. The same exercise for i x = 42° results 
in a much better overall agreement between the data and sim- 
ulation, with the simulated radial points matching the rough 
positions and scatter of the real ones. At i x = 50°, the sim- 
ulated views to the north (corresponding to the data viewed 
close to 43° and 37°) again fail to correctly place the radial 
material. Specifically, the simulated radial material viewed 
from the north at 6.5° looks more like echo data at i = 37°, 
whereas it should match up with the view at i = 43°. Within 
its limited predictive capacity, these results favor a true incli- 
nation in the vicinity of i x = 42°. 

The model can also test whether fitting ellipses to these par- 
ticular data can provide meaningful results. For each panel in 



Figure|25] we fit ellipses through all points, and through only 
the simulated shell points. Those fits, as well as their elliptic- 
ities (bja), are indicated in each box. The actual ellipticities, 
from left to right, should be 0.96, 0.99, 1.0, 0.99, and 0.96, 
with the centroid's y-offset changing roughly as shown in the 
top row. 

In very few cases does either fit give the correct ellipticity, 
and only in one-third of the cases are the offsets correct. Most 
fits favor an orientation with the major axis along y' . Con- 
sidering only the face-on view (middle column) the fits have 
an error of 3-5% including all data, and 3-14% using only 
the shell points. In the top two rows, the fits are, on aver- 
age, good to 5%, with considerable room for additional error. 
The bottom panel shows highly-eccentric fits with large cen- 
troid offsets yo. These are "false positives," since x 2 nas been 
minimized, but the ellipse does not reproduce the actual cross 
section. Note however that in all panels, the centers rarely 
deviate east or west, thus a measured offset in Xo may be gen- 
uine. 

Since the uncertainty in a fit is much larger than the change 
in ellipticity with inclination, it is difficult to perform an anal- 
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FIG. 26. — Similar to Fig. 1191 but for the favored orientation of i x = 40° 
and i y = 8°, and showing \z'\ along the abscissa. Overplotted is the average 
radial profile, measured as the average position (and standard deviation) of 
points binned along |z'|. 



ysis similar to that in 36.1.3l using ellipse fitting. Nonetheless, 
we fit a variety of ellipses to the cross-sectional profiles of 
the CD echoes for i x = 40° -42°. We consistently find the 
cross section's center is shifted east by ~ 1.5 ly. This i s con- 
sistent with the offset in shell material first noted in 36.1.11 
At these inclinations, cross sections also have an ellipticity of 
b/a = 0.95 ± 0.05, with the major axis aligned toward north. 
However, as noted above, this is not a robust measurement for 
this particular dataset. 

6.2. A Complete Picture of the CD 

The analyses of the previous subsections suggest the CD 
echoes lie on a surface of rotation about a single axis, with 
a favored orientation of i x = 40° and i y = 8°. The radial pro- 
file for this orientation is shown in Figure[26] All features are 
quantitatively very similar to those in Figure^] Shell points 
a re desc ribed by nearly the same linear fits as for i x = 37° 
( 36.1.2> . while the radial and jet points are fit by slightly- 
steeper lines. Both the radial profiles and ellipse-fitting sug- 
gest that this structure's axis may be offset to the east of the 
SN by < 1 ly. This is not a robust result, and since such a 
small offset will have a negligible effect on the shape of the 
CD structure, we ignore it for now. 

Assuming these data also trace a structure that is symmetric 
about z' = 0, an average radial profile of (it 1 ) versus \z' | can be 
constructed by binning the points along the axis. The average 
radial profile (and standard deviations in each bin) is shown 
over the data in Figure [26] 

The probable CD structure is visualized by revolving this 
function about the symmetry axis, reflecting it about the equa- 
tor, and reinclining it to the favored orientation, as shown in 
Figure I43h-l>. Panels (c)-(d) show this complete structure 
plotted in monotone grey, overlaid with the actual echo points 
from Figure [21 For these panels, shell points are colored 
blue, while radial points are red and the jets are green. 

This structure is fairly complicated, and can not be de- 
scribed by a simple geometric function. Figure \^\ shows a 
scaled cartoon of the salient features. The "shell" echoes ap- 
pear to lie along a prolate spheroid with a polar axis of 20 
ly and equatorial axes of 11 ly. However, the ends of this 
spheroid have been extended into tapering cones with open- 
ing angles of about 35° (the "radial" echoes), extending from 
16 to 28 ly from the SN. Embedded within this prolate struc- 
ture are narrower, tapering cones (the "jet echoes") extending 
from 10-26 ly from the SN, with an opening angle of about 




FIG. 27. — C artoo n sketch of the salient structures traced out by the CD 
model from Fig. 1261 This figure is to scale, with the length of the two orien- 
tation arrows equal to 10 ly. Solid lines trace the average radial profiles, and 
the scatter above those profiles are indicated by the solid greyscale regions 
surrounding each line. 



20°, and a maximum radius of 3 ly. Some shell points in 
Figure |26] are located at n 1 ~ 15 ly, and appear disconnected 
from the rest of the material. These points, which we call the 
"spurs," lie along a cylindrical annulus that smoothly encir- 
cles the CD. 

It is unclear whether the prolate shell is continuous in z' 
along its equator (dashed equatorial line in Figure l27l . since 
no echoes were observed from that region. Whether this is due 
to shadowing from material closer to the SN is addressed in 
^9] The northwestern- WFPC2 echoes are positioned at z' ~ 4 
ly along the axis, and at roughly 8 ly in radius. This could 
suggest the prolate shell is pinched at its waist to a smaller ra- 
dius of 8 ly. However given the very limited spatial sampling 
of these inner echoes, there is little evidence that they lie on a 
uniform structure. 

At this point in the analysis, evidence that the structure 
containing these echoes has rotational and mirror symmetries 
is indirect, and since construction of the probable CD neb- 
ula shown in Figure|^]relies on these two assumptions, they 
merit a brief discussion. First, we recall that Figures [T8l and 
1 191 are consistent with surfaces of rotational symmetry. How- 
ever, rotational symmetry need not imply mirror symmetry, as 
in the cases of a paraboloid or cone. As we will see in the next 
two sections, light echoes from material interior to the CD 
nebula lie on axisymmetric surfaces which do exhibit mirror 
symmetry. We therefore make the (reasonable) assumption 
that stellar winds producing symmetric structures at small dis- 
tances would also produce such structures at larger ones. 

In con cluding this section , we first comment on the sugges- 
tions in 36. 1.31 and 36. 1.41 that the CD has an inclination of 
i x = 50°. This is a result of the particular portion of the CD 
that was imaged in echoes, rather than an intrinsic property of 
the CD itself. The probable CD is rendered in Figure|28] now 
shade d by the relative brightness of dust scattering ( Sugerman 
120031) . wich assuming uniform density and composition is a 
function of distance and scattering angle. Panel (a) shows 
what the entire CD would look like though a complete set of 
observations well-sampled in time. Panel (b) shows only that 
portion of the complete CD that could have been observed in 
our epochs of observation. When viewed from the CD sym- 
metry axis (panel c), we see that brighter echoes form a sharp, 
circular ridge that could easily be mistaken for delineating the 
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FIG. 28. — Simulated CD showing relative brightness of echoes, where 
greyscale indicates relative brightness, with darker grey corresponding to 
brighter echoes. To distinguish faint points, they have been circumscribed 
with darker grey circles; the reader is reminded that they shade of grey inside 
these circles indicates brightness, (a) Full CD, viewed face-on. (b) Random 
points illuminated by the echo parabolae from all epochs, (c) Panel (b) in- 
clined along the CD axis of symmetry. 
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FIG. 29. — Spherical positions of NH echoes (see Fig. 1181 . NH-norfh are 
plotted with open circles, and NH-south with filled circles. The best-fit lines 
through the horns (solid line) and the other data (dotted line) are indicated. 



sharp edge of a narrow cylinder that was seen from i x = 50° in 
Figure l20l c). 

Finally, if the CD is indeed mirror symmetric, why were 
echoes not seen from the far half of the nebula (material at 
— Z 7 )? In particular, there is a small portion of material north 
of the SN interior to the late-time paraboloids (northern points 
jus t interior to the large parabola in Figure fl5? fc)). Figure 
\28l b) shows that echoes of this small northern portion were 
very faint due to the large scattering angle, and hence were 
unlikely to be detected. On the other hand, a segment of the 
spurs along the closer half of the CD should have created very 
bright echoes at a radius of p = 7 ly, roughly from PA 315°- 
30°. Very little signal was detected in this region, which may 
suggest that the spurs are isolated clumps, rather than a uni- 
form feature. 

7. NAPOLEON'S HAT 

In contrast to the previous case of the CD echoes, there is 
much less ambiguity over the possible structures of which the 
NH echoes are a part. As shown in Figurefl5t/-/'. NH appears 
to be some surface of revolution with a well-defined axis, such 
as a prolate ellipsoid, cylinder or frustum. However, it is not 
immediately clear whether the two horns are part of such a 
structure. If not, where do they lie, and if so, why has such 
a limited portion been illuminated? The geometry of the NH 
echoes will be quantified in the following subsections. 

7.1. Geometry ofNH Echoes 
7.1.1. Spherical Coordinates 




1*1 (iy) 

FIG. 30. — Cylindrical positions of NH echoes (see Fig. 1191 for th e bes t-fit 
inclination of 40° south and 6° 5 east. Symbols are the same as Fig. [29] We 
assume NH is symmetric about the equator, and plot z'|. The solid line is the 
best-fit through all points. 



As with the CD, we examine the NH echoes in spherical 
coordinates, shown in Figure[29] The bulk of the echoes fill a 
horizontal region in <f>, with points distributed between n' = 3 - 
6 ly, but concentrated around r = 4-5 ly. In these coordinates, 
a thick spherical shell with radius r ~ 4.5 ly and width Ar = 1 
ly would populate the same region. However, such a structure 
is unlikely: if the data lay along a spherical shell, echoes in 
front and south (or behind and north) of the SN should have 
been observed. 

The horns are also inconsistent with a spherical structure, 
as they are instead located in a highly inclined, almost ra- 
dial region extending from 6-1 1 ly from the SN. Rather, the 
total ensemble of NH echoes can be understood to lie along 
a cylinder or hourglass. The distributions of points on such 
structures, viewed in spherical coordinates, populate a hori- 
zontal region in <p near the equator, while points far from the 
equator will be highly inclined. 

7.1.2. Cylindrical Coordinates 

As with the CD echoes, studying points in cylindrical co- 
ordinates requires knowing the axis z' ■ Since the candidate 
structures for the NH echoes are an hourglass or cylinder (a 
specialized case of an hourglass), we start by fiting a biconi- 
cal f rustum to all NH echo points, as explained in Appendix 
IB. 41 We first restrict the figure's cross section to be circular, 
and the SN center to lie on the symmetry axis. With these 
constraints, the best-fit structure (RMS scatter of 0.73 ly) is 
an hourglass with waist radius of 3.7 ly and opening angle 
36-39°, inclined i y ~ 40° south and i x = 2-4° east of the line 
of sight. Allowing for an elliptical cross section improves the 
X 2 by nearly a factor of two (RMS=0.55 ly), yielding semi- 
major and minor axes of 4.6 and 3.8 ly, with opening angles 
at the axes of 28 and 40°, respectively. The major axis is ro- 
tated 96° ± 4° east of north, and the axis of the hourglass is 
oriented with i y = 39.°5 ± 0.°5 south and i x = 6.°5 ± 0.°5 east of 
the line of sight. This fit is robust to perturbations in the fitted 
parameters, including spatial offsets of the waist from the SN. 

Figure|30|shows the radial profile of these points as a func- 
tion of absolute position along the axis (this assumes the data 
have reflectional symmetry about z' = 0. This figure is con- 
sistent with an hourglass shape. The best fit through all data 
is given by tt' = (4.13±0.06) + (0.42±0.04)|z'|, which cor- 
responds to a waist radius of 4. 1 ly and an opening angle of 
46°. While there is a large scatter (RMS=0.5 ly) about this 
line, note that there are similarly-sloped loci of points < 1 ly 
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FIG. 31. — Cross-sectional view of NH echoe s lo oking down the best-fit 
inclination axis. Symbols are the same as Fig. 1291 The best-fit ellipse is 
indicated by the solid line. 



above and below it. This is further suggestive of a structure 
with elliptical cross-section. 

The horns extend along the axis from 1-7 ly in front of the 
SN, and radially 4-8 ly from the axis. They are best-fit by 
the line tt' = (4.4 ± 0.2) + (0.52 ± 0.07)z', which corresponds 
to a frustum with inner radius 4.4 ly and an opening angle 
of 55°. The horns do not appear to be truly linear, but curve 
to shallower slope at larger z', as if they lie on a surface of 
parabolic, rather than linear, profile. 

We directly examine the ellipticity of the NH cross sec- 
tion by obbserving the distribution of points looking down 
the be st-fit a xis of inclination (40° south, 6?5 east). Simi- 
lar to 36. 1.41 with the CD, echo positions (x',y') in the plane 
orthogonal to z' (see Figure II 21 are normalized by the aver- 
age radius tt' as a function of z' given by the best-fit line in 
Figure [30] The resulting distribution, shown in Figure I3T1 is 
fit with a 5-parameter ellipse, as described in Appendix IB. 21 
The best-fit ellipse has negligible offset from the origin, and 
b/a = 0.82 ±0.02 with the major axis rotated 103° + /-3 east 
of north. This is consistent with the hourglass fit above, which 
yielded b/a = 0.83 and major-axis rotations of 92-101°. 

We use these values to remove the ellipticity from the ra- 
dial profile, yielding the distribution shown in Figure [32] As 
expected, the horns have shifted to smaller radii, and the scat- 
ter between the rest of the NH points has diminished. Also 
shown is the average radial profile of the NH echoes, gen- 
erated by binning the points along the inclination axis and 
computing the average radius about that value. The implied 
structure appears cylindrical with a semi-minor axis of 4 ly up 
to \z' I = 1 .5 ly, after which it has a fairly constant slope. The 
best-fit line through the entire profile (solid line) is given by 
tt' = (3.53 ±0.23) + (0.52 ±0.06)z', corresponding to an hour- 
glass with semi-major and minor waist axes of 4.4 and 3.5 ly, 
and an opening angle near 54°. Excluding the horns, the fit 
is considerably shallower (dotted line), suggesting an opening 
angle of 36°, and axes of 4.8 and 3.8 ly. Despite the ambiguity 
in functional form, these fits are all consistent with the values 
found by the hourglass fitting. 

7.1.3. What are the horns? 

If the horns are part of the same hourglass-like structure 
containing the rest of the NH echoes, why was such a limited 
section of the much larger surface illuminated? Assume the 
horns lie along the solid line in Figure|32] This line can be re- 
volved about the z' axis, then reinclined from the line of sight 
by the best-fit NH orientation to generate a complete frustum, 
as shown in Figure l33b. As in Figure [28] grey scale indicates 
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FIG. 32. — Data from Fig. 1301 with the ellipticity from Fig. l31l removed. 
Overplotted is the average radial profile, measured in bins along \z'\. Single 
points have fixed errors. The best-fit line through the complete profile is 
indicated by a solid line, and the best-fit through points excluding the horns 
(|z'| < 3 ly) by the dotted line. 
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FIG. 33. — As Fig. l28l but showing the frustum containing the horns. Darker 
shading indicates increased surface brightness of an echo, (a) The complete 
frustum as revealed through a complete, regularly sampled set of echoes. 
Random points along the intersection of the frustum with echo parabolae 
from our observations are shown in panels (b) (face-on) and (c) (from the 
top). Greyscale is normalized to the brightest and faintest points plotted, thus 
shadings differ between panels (a) and (6)-(c). 



relative brightness of echoes based on their scattering angle 
and distance from the SN. Points along the northern limb are 
the brightest, since the scattering angle is lowest. 

Random points on this frustum along echo paraboloids ob- 
served in our data are shown in panels (b)-(c). The northern 
limb was never illuminated by echoes as it lies inside the ear- 
liest paraboloid. The earliest intersection occured just south 
of this limb, and over time, increasingly-southern slices were 
illuminated, becoming fainter with increasing epoch. The 
brightest echoes came from the first few epochs, appearing 
as linear, radial features at roughly constant north-south posi- 
tion. Furthermore, at the earliest epochs, the illuminated ma- 
terial is oriented close to the line of sight, thus these echoes 
appeared even brighter due to a limb-brightening effect. Thus 
the limited portion of this structure that was illuminated as 
"the horns" is consistent with these points lying on the same 
hourglass as the rest of the NH material. 

7.2. A Complete Picture of the NH Structure 

The probable structure containing the NH echoes is an 
hourglass, the waist of which has an ellipticity of b/a = 0.82, 
with a semi-major axis of ~ 4.9 ly rotated about 102° east of 
north. Bootstrap analyses of these fits yield errors of about 3% 
for each parameter. If the horns lie on a surface with similar 
eccentricity, the hourglass flares to a semi-minor axis of ~ 7 
ly about 7 ly from the SN along the axis of symmetry. The 
entire structure is oriented at i y = 40° south and i x = 7° east 
of the line of sight. To visualize this, the average radial pro- 
file from Figure[32]is reflected about the equator and revolved 
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FIG. 34. — As Fig. 1271 but for the NH echoes, with the orienta tion arrows 
measuring 4 ly. The solid line traces the radial profile from Figure l32l and the 
width of the greyscale region traces the scatter of points about that profile. 



around the axis, then reinclined from the line of sight by the 
inclinations just given. The resulting distribution is rendered 
in Figure \43e-h. 

A much larger portion of this nebula was sampled by echoes 
than was actually observed. This was discussed in some detail 
above within the context of the horns, from which we con- 
cluded that the majority of the signal from the near half of 
the hourglass was very faint due to the larger distance and 
scattering angle of that material compared to the dust at the 
waist. The same argument applies to the far half of the hour- 
glass, explaining why no discrete signal was observed from 
the northern parts of the structure at later times. 

A number of small-scale features not represented in this fig- 
ure are better discussed in the cartoon sketch of the NH echoes 
shown in Figure [34] where grey-shading represents the scat- 
ter in points about the average profile. In the first two bins 
(|z'| < 0.5 ly) in Figure|32] material extends beyond the aver- 
age scatter to smaller and larger radii of tt' < 3 and tt' > 4.5 
ly. Furthermore, over one-third of the NH points lie within 
these first bins, which implies these data trace out an equato- 
rial enhancement in both density and radial position. This is 
shown in Figure [34] as small spurs pointing toward and away 
from the SN in the equatorial plane. A similar, but less sig- 
nificant, inward-pointing feature is seen roughly 1.5 ly along 
the axis. CKH95 reported a small set of echoes within 3 ly of 
the SN that are coincident with the ER plane, which they in- 
terpreted as evidence for an extended circumstellar equatorial 
disk. The NH equatorial echoes appear to extend this disk to 
a radius of 5 ly, with a thickness of > 0.75 ly. 

7.3. Comparison with Previous Work 

As NH was the subject of interest in the early 1990's, we 
compare the newly revealed m aterial with previous interpre- 
tations. Wamp ler et al.l dl990al) first reported the discovery 
of the "Napoleon's Hat" echo in narrow-band [O III] and 
[N II] imaging as diffuse nebulosity, which blended smoothly 
into background filamentary structure that defines the east- 
ern and western ri ms of a "V"-shaped "d ark bay" south of 
the SN (see Fig.[6]i. Wam pler etaD lll990bl) further noted that 
the larger, arcm in-radius echoes from interstellar dust (e.g., 
IXu et alj [l995) fade when crossing over this bay. A natural 
interpretation was that this bay is an evactuated region, and 



no light echoes were seen there because of its particularly low 
density. As the background filaments loop around the SN, 
they proposed that the progenitor's MS winds evacuated the 
bay, and its V-shape can be attributed to the proper motion of 
t he progenitor through the ISM. 

Wang & Wampler ( 1992) reanalyzed their ground-based 
data and were able to trace the NH echoes as early as day 
890, but the echoes had faded by day 1650. An interpre- 
tion that these echoes lie along a disk coplanar to the inner 
ER was disfavored, since the flux had the wrong distribu- 
tion and shape for a solid, dusty disk, and since the echoes 
mo ved very little in time. Rather , Wang & Wampler ( 1992) 
and Wang. Dvs on. & KahnlJT9 93) considered these structures 
consistent with a parabolic bow-shock to the north of the SN, 
which could be further interpreted as the terminal edge of the 
dark bay evactuated by the MS progenitor. 

Given the limited data at hand, and considerable confusion 
from the background filaments in the LMC, the bow-shock 
model was an attractive scenario. However, the current anal- 
ysis has revealed a much more extensive volume of material 
illuminated by the SN, the 3-D positions of which are incon- 
sistent with a bow shock. Rath er, these data imply that the dis- 
favored interpretation of Wang & Wampler ( 1992)-that NH is 
in the equatorial plane-was correct. 

Suppose the star does have a large proper motion, and that 
the dark bay is the region evacuated by the MS wind. What- 
ever the star's proper motion through the ISM, it must have 
moved at least 25 ly north during its lifetime to create this bay 
and end at its observed position. The CD echoes from this re- 
gion have illuminated a prolate structure that ends at roughly 
the same position on the sky as the center of this bay. This 
is also inconsistent with the bow-shock model. Most likely, 
the filamentary nebula is not associated with the SN, and we 
propose that the correlation between the position of the NH 
echoes and the "V"-shaped nebula is a positional coincidence. 

8. CIRCUMSTELLAR ECHOES 

As noted in §l\ the inner 3" surrounding SN 1987A have 
been previously studied by CKH95 using a subset of the data 
presented here. Employing PSF-subtraction and rudimentary 
PSF-matching, the authors detected echoes along a bipolar 
nebula with an hourglass shape, in the sense that a narrowed 
waist separates two open-ended lobes. A smaller set of echoes 
coplanar with the ER were identified as evidence of an ex- 
tended circumstellar disk. 

The current analysis differs from CKH95 in a number of 
ways. In the previous work, the authors payed particular at- 
tention to the inner few arcsec around the SN, and their re- 
duction may have less noise from PSF-subtractions of the 
SN and s tars 2 and 3 than from our automated reduction 
pipeline ( 34.3.3> . The general PSF-matching algorithm has 
been greatly improved since 1995, and we have implemented 
a number of custom modifications to minimize difference- 
residuals and to increase the detection threshold of echo sig- 
nal. Since we were interested in echoes from a much larger 
range of radii, we also have increased the spatial extent over 
which circumstellar (CS) echoes have been sought. Finally, 
we used a very different technique for identifying the posi- 
tions of echoes. Sufficient differences exist between the two 
methods that the current data merit an independent analysis. 

8.1. Geometry of CS Echoes 

An ex plora tion of the geometry of these echoes (summa- 
rized in 38. 2> is facilitated by the work of CKH95, which es- 
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FlG. 35. — Cylindrical profile of tt' versus \z'\ for the CS ec hoes inclined 
by i x = 45° and i y = 8°. Symbols correspond to points in Fig. 1 15b /. with 
closed and open circles denoting cyan and yellow extended flux points, and X 
marking the red hourglass. Points below the solid line are potential members 
of the CS hourglass, and comprise d istribu tion "A," while those also under 
the dashed line are distribution "B" ( 88.1.11 . 





FIG. 36. — Same as Fig. l43fr -7i, but showing only those CS p oints believed 
to lie along the hourglass (the case B distribution from Fig. 1351 . 



TABLE 2 

Observed Geometric Parameters of the 
ER 



tablished the cylindrical symmetry axis of the double-lobed 
CS hourglass. We cannot immedi ately e xplore such a fea- 
ture in our data since, as noted in 35.4.31 it is unclear which 
extended-flux echoes are associated with this hourglass. How- 
ever, CKH95 reported a small number of echoes along this 
southern lobe, and we can use their data to help disentangle 
our points. 



Parameter 


Unit 


Value 


Centroid Offset 


(mas W) 


18.7 ±0.4 


Centroid Offset 


(mas S) 


2.0±0.7 


Inclination Angle 


(°) 


43.0 ±0.1 


Semi-Major Axis 


(") 


0.820 ±0.005 


Semi-Minor Axis 


(") 


0.600 ±0.005 


PA. of Major Axis 


(° E of N) 


81.1 ±0.8 



8.1.1. Geometry in Cylindrical Coordinates 

We begin by fitting a biconical hourglass to the C KH95 
data, using the \ 2 minimization described in Appendix IB. 41 
Assuming the frusta have circular cross-secitons, the data are 
best-fit by a structure inclined by i x = 45° south and i y = 8° 
east of the line of sight. This compares favorably to the ori- 
entation reported in CKH95 using a minimization-of-scatter 
estimator, but our rotation i y is larger than theirs (3° ± 5°). 
We adopt this as the preliminary orientation axis, and trans- 
form all echo points into cylindrical coordinates (see Figure 
\12i about this axis. Assuming the structures are symmetric 
about the equator (zf = 0), we plot in Figure[35]the radial pro- 
file of tt' as a function of absolute distance \z'\ along the axis. 
The region containing northern-hourglass echoes (marked by 
x ) includes some extended-flux material, which may corre- 
spond to additional points along the northern and southern 
lobes. 

The northern-hourglass echoes in Figure[35]can be consid- 
ered as lying along one of two regions. Case "A" is roughly- 
delimited as points below the solid line, which follows the 
outer edge of the echoes at any given point along the axis. 
When binned, the northern-hourglass points follow a some- 
what curved trajectory that approaches constant radius with 
increasing \z'\- Such a region is approximately bounded by 
the dashed horizontal line, which we denote case "B." While 
a few extended-flux points straddle the solid boundary, the 
majority are located above it, suggesting that points above the 
solid line are part of a different structure, and only points be- 
low it can be part of the hourglass. We proceed by testing 
which distribution of points is most-likely on such a structure. 

Biconical hourglasses are fit to the case A and B distri- 
butions. Case A points are consistent with an hourglass in- 
clined 41°-42° south and 10°-13° east of the line of sight. 
The fit to distribution B is similar, with an inclination of 40°- 
41° south, 12°-15° east. This fit has a \ 2 nearly three times 



smaller than for case A, while there are only 16 fewer points, 
suggesting that the points above the dashed line in Figure [35] 
are not directly part of the hourglass structure. Accordingly, 
we choose the B distribution as the probable set of hourglass 
points. These points are rendered in Figure [36] As expexted 
from the earliest echo parabolae, only a small arc along the 
southernmost limb of of the southern lobe is present. 

The best-fit hourglass to these data is inclined at (i x ,i y ) 
of (-40.7° ±0.2°, 14.7° ±0.7°), with a westward shift of 
0.04 ± 0.01 ly. The waist has an ellipticity b/a = 0.94 ± 0.03 
and semi-major axis 1.01 ±0.2 ly, rotated 10° ±1° north of 
east. As noted in Appendix 151 these errors are the scatter in 
parameters fit to random subsets of the data. 

8.1.2. Constraints from the Equatorial Ring 

If we assume the three-ring nebula lies symmetrically about 
the hourglass, then the positions of the rings on the sky should 
intersect, or be symmetrically bounded by, the hourglass. 
These rings are commonly assumed to be planar circles, ap- 
pearing elliptical due to their inclination. As such, their ob- 
served positions should be reproduced by simple geometric 
figures viewed in the hourglass frame of reference. Unlike 
the ORs, for which only limited velocity data exists, the ve- 
locity shear across the ER indic ates it is a planar structure 
(Crotts & Heathcote 1991, 120001) . and its projected geometry 
is very-well constrained (see below). The ER is therefore the 
best candidate for t his ex ercise. 

Sugerman et al. (2002) measured the ER geometry by fit- 
ting ellipses to a high-resolution WFPC2 image. The reader 
is referred to that paper for a thorough description of the 
method. Briefly, the ER was divided into bins of equal ar- 
clength, and the flux in each bin was normalized to the sec- 
ond brightest pixel value. Ellipses were then fit to the bright- 
est N% of the pixels in each bin, where N varied from 50% 
to 95%, and the resulting parameters averaged to give final 
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FIG. 3 7 . — Drizzle-recons t ructed (0."23 pix" 1 ) WFPC2 F656N image of 
the ER from Sugerman et al. i2002). The best fit ellipse to the observed ER 
(Tablel2l is indicated with the thick dotted line. The projected position of an 
ellipse oriented by the PYR triplet (-16?5,-14°,41?5) is shown with a thin 
dashed line, corresponding to the best fit to the hourglass from 88.1.11 The 
projected view of the ellipse that best-matches the ER from 88.1.2l oriented by 
(-9°,-8°,41°) is indicated by the thin solid line, which nearly exactly matches 
the thick dotted line. 



values. However, their choice of fit ting algorithm was not as 
robust as that described in Appendix lB.21 We reproduced their 
fit by substiting our improved minimization algorithm, the re- 
sults of which are given in Tableland drawn with heavy dots 
in Figure l37l 

If the ER is a planar ellipse, then viewed at the hourglass 
orientation, it should overlap closely with our best-fit. For 
all the orientations listed in the previous section, no combina- 
tion of ellipse parameters (offsets, axes and major-axis rota- 
tion) can fully reproduce the observed geometry. The best fit 
achieved using those orientations is shown with a thin dashed 
line, and has an average deviation of 18 mas. This error is 
more than double the uncertainty of the ER axes, and is as 
large as its centroid offset, thus the orientation provided by 
the hourglass fit is not ideal. 

Allowing the orientation of the ellipse plane to vary as well, 
the absolute best fit to the observed ER deviates by only 3 
mas. Shown by the thin, solid line, this ellipse has an orienta- 
tion of i x = 41° south and i y = 8° east, and with the major axis 
rotated 9° north of east. The ellipse has a major axis 0"82, 
bja = 0.98, and the centroid is shifted 19 mas west of the SN. 
A circle at this same orientation fits the ER with three-times 
larger % 2 . Provided the ER is connected to the CS hourglass, 
this is the first direct measurement of the ER's orientation and 
deprojected geometry. 

8.1.3. Geometry in Cylindrical Coordinates II 

Why do the hourglass and ER fits yield values of i y dis- 
crepant by over 6° ? The hourglass is well-sampled only along 
the southern edge of the northern lobe, otherwise, the data 
are sparsely sampled. This distribution is sufficient to con- 
strain the north-south inclination i x , but allows great freedom 
in east-west rotations. We test the newly-proposed orientation 
by repeating the hourglass fit, this time holding i y fixed at 8°. 

The result is an hourglass inclined i x = 41.°2±0.°2 and 
shifted 0.1 ±0.02 ly west; a waist with semi-major axis 1.04± 
0.02 ly rotated 9° ± 1° north of east ; and bja = 0.94 ±0.01. 
The same caveat for the errors from IHTUappl ies. The frusta 
have a half-opening angle of 12° ± 1, which is close to that 
found by CKH95. This is very similar to the orientation sug- 
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FIG. 38. — Same as Fig. 1351 but with data inclined by the favored CS 
orientation of i x = 41° and i y = 8° and an ellipticity of 0.94 removed from 
the radii. Circles that are now considered part o f the hourglass are marked 
by open squares, and the case-A only points (Fig. 135 1 are marked by pluses. 
The average radial profiles to each structure are o verplotted. The dashed line 
shows the inner edge of the belt, used in 418.31 Note that there are many 
binned averages in this plot, which are discussed in the 88.21 



gested by the ER fit, and therefore we adopt it for the CS ma- 
terial. Note tha t fitting an hou rglass confirms the findings in 
SI8.1.2l and by Sugerman et al. (2002) that this nebula is offset 
slightly west of the SN. 

The cylindrical-radial profile for the data oriented about this 
new axis is shown in Figure[38] An ellipticity of bja = 0.94 
has also been removed from each radius, based on its posi- 
tion angle along the elliptical hourglass fit given above. Ex- 
tended flux points that are now associated with the hourglass 
are plotted as squares, and the points belonging only to the 
case A distribution are shown with plus marks. A large subset 
of the extended-flux echoes (circles) are concentrated within 
1 ly of the equatorial plane (z' = 0), between 1 .5-4 ly from 
the SN, and appear to form a thick equatorial waist or belt cir- 
cumscribing the hourglass. As in the previous two sections, 
we measure the average radial profile of these points in bins 
along \z'\- There is much structure to these data, which re- 
quires the fitting of six profiles. These, and their interpreta- 
tion, are discussed in the next section. 

8.2. A complete picture of the CS nebula 

Surrounding SN 1987 A is a complex CS nebula, inclined at 
i x = 41° south and i y = 8° east of the line of sight, translated 
< 0.1 ly to the west of the SN, with an elliptical cross-section 
(bja~ 0.94) whose semi-major axis is rotated 9° clockwise 
from east. We visualize this nebula by reflecting each radial 
profile about the equator, revolving it around the cylindrical 
axis, then reinclining it by the best-fit orientation. The ellip- 
ticity removed in making Figure [38] was reintroduced when 
revolving the profiles about the cylindrical axis. The result, 
rendered in Figure E3j-j, shows the probable structure of the 
CS material. Colors are discussed below. The next two panels 
(k-J) compare the complete structure in grey to the observed 
echoes. Here, colors are similar to Figure lTBp -/. only now all 
hourglass points are colored red, and the Case-A only points 
(Fig. [35} are colored green. For completeness, we show only 
the CS hourglass in panels (m)-(p). 

Figure [39] shows a scaled cartoon of the relevant CS fea- 
tures. In the following subsections, we discuss each feature 
based on its radial profile and three-dimensional shape. 
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FIG. 39. — As Fig. 1271 but for the CS echoes, with the orientation arrows 
measuring 2 ly. The expected positions of the ORs are indicated by thin 
dotted lines f flO . 



The innermost radial structure is the hourglass, shaded red 
in Figurel43l-/ and shown by itself in panels (m)-(n). In Fig- 
ure |38j these are the crosses (x) and squares within about 
tt' < 1 .5ly. This is fairly cylindrical, with a semi-minor axis of 
1.2 ly, tapering close to the equator and opening just slightly 
at large zf- The profile is well-described by a line of slope 
Air' I Az' = 0.14 and zeropoint tt' = 1 .02, suggesting an open- 
ing angle of the wall of only 8°. A very small number of 
points suggest that it also tapers around z! = 1 .5 - 2 ly, but this 
feature is marginal at best. 

A number of small-scale features are worth noting in the 
hourglass. The northern-lobe points extend to a relative max- 
imum in radius (tt' ~ 1.5 ly) at = 1.4 ly, tapering slightly 
in radius on either side; a similar feature is also seen around 
\z'\ = 0.9 ly. In Figure[38] many of the extended-flux points lie 
just outside the average position of the hourglass. Since most 
of these points are part of the southern lobe, this could sug- 
gest the southern lobe is wider than its northern counterpart. 
A more likely interpretation is that many of the southern-lobe 
points represent dust just outside the densest part of the hour- 
glass, which was illuminated prior to our earliest epoch (see 
Fig.Efc). 

Consistent with CKH95, we do not see evidence for a "cap- 
ping surface" to this structure, as would be expected from 
a bipolar peanut-like nebula, the result of many interacting- 
winds formation models. To distinguish this hourglass from 
the larger one containing the NH echoes, we refer to each by 
their association-i.e., the latter is the "NH hourglass." 

As in previous sections, we address why only a limited 
amount of this complete hourglass structure was illuminated 
by echoes. Referring to Figure l43lp ). only the southernmost 
portion of the near (or south) half of the hourglass was probed 
by the earliest echo. While the greater majority of the far (or 
nothern) half of the CS hourglass lies inside the earliest and 
latest echo parabolae, the particular sampling of our obser- 
vations probed the northern edge only at the earliest epochs. 
A comparison of this panel with Figure l36T £>) shows that we 
have detected roughly all the material that could have been 
observed. 

8.2.2. The Equatorial Plane and Belt 

Immediately surrounding the CS hourglass is the thick 
waist of extended high-surface brightness, which we call the 
"belt," and have colored green in Figure 03}-/ It extends ±1 
ly along the zf axis, and from 1.5 to 2.5 ly in radius. To better 



represent the extent of this material, two profiles have been fit 
in Figure[38] binned both in radius and along the axis. Beyond 
the outer radius of this belt, echoes lie fairly well constrained 
to the equator, tracing a thinner (0.5 ly thickness in z') but ex- 
tended equatorial plane to an outer radius of tt' = 4 ly. This 
feature has been colored blue. The belt does not appear to ta- 
per smoothly into the thinner equatorial plane, but rather has 
a fairly sharp transition at its outer boundary. 

8.2.3. The Inner and Outer Walls 

The rest of the CS material lies along one of two "walls." In 
Figure|38] the innermost of these is marked mostly by plusses 
at \zf\ > 1 .5 ly, while the outer wall is marked with circles 
from 0.9 < \z'\ < 1.9 ly. The inner wall extends 1.5 to 2.6 ly 
along the axis, and may be considered to join with the hour- 
glass around \z'\ = 1.5 ly. The profile is positively sloped, 
which translates into a wide cone with half-opening angle of 
40°. This may simply trace a decreasing dust-density profile 
extending out from the sharp peak of the hourglass. A single 
datum is located at \z'\ = 3.4 ly, and it is unclear whether this 
can be trusted as a real feature. 

The outer wall bridges the belt and the inner wall at radii 
between 2.6 and 3.4^1.0 ly. These have been colored cyan 
and gold in Figure 143 M. Unlike the inner wall, the average 
profile is concave, reaching its maximum radius roughly 1 .4 
ly along the axis. This corresponds to a local maximum seen 
along the hourglass, and to the expected positions of the ORs 
iCrotts & Heathcotel bOOO). Noting the conspicuous lack of 
echo signal between the hourglass and outer wall at this posi- 
tion, we question whether this is indicative of episodic mass 
loss, and if the OR formation mechanism is somehow related. 

8.3. The Geometry of the Outer Rings 

Whil e great attention has been placed on the geometry of 
the E R (Uakobsen et al.lll991tlpTaitet^fT9^5l |BuiT ows et al. 
119951) . relatively little work exists on the ORs. In 38. 1.21 the 
ER was used as a constraint to the geometry of the CS nebula. 
Now that a three-dimensional map of the nebula exists, this 
problem can be inverted to ask instead how the echoes can 
c onstrain the 3-D geo metry of the ORs. 

Burrows et al. (1995) reported the NOR is well fit by an 
ellipse with b/a = 0.678, and a semi-major axis of 1 ."77 ori- 
ented at RA. 70.°7, while the SOR is fit by an ellipse with 
b/a = 0.513, and semi-major axis 1 ."84 oriented at RA. 90. 
Assuming both rings to be circular in their reference frames, 
they are inclined roughly 43° and 31° to the line of sight, re- 
spectively. Accordingly, if the SOR has the same inclination 
as the ER and NOR, it cannot have a circular shape. They 
also note that the segment connecting the centers o f both ORs 
misses the SN and ER centroid by ~ 0."4. ICrotts & Heathcotd 
(2000) use these orientations to fit a kinematic model to the 
ring velocities as measured in STIS and ground-based echelle 
spectra. This model suggests the NOR and SOR are located 
1.32 and 1.04-1.24 ly (depending on whether one forces the 
SOR to be circular or not) from the SN along the inclination 
axis. 

Figure |40| shows the complete three-ring nebula from the 
same WFPC2 image used to study the ER in 38.1.21 We mea- 
sured the positions of each OR at various RA.s by perform- 
ing the same radial-profile fitting as that used in measuring 
echoes. With the positions of the rings on the sky known 
to better than one pixel (0."023), the only unknown for each 
point is its line-of-sight distance from the SN. If one assumes 
the ring lies on or about the CS hourglass, each z position can 
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FIG. 40.— The full WFPC2 image of the three rings, from which Fig.BTI 
is taken. To increase the display range, brighter pixels are displayed with 
a separate color stretch. Heavy crosses show the positions of the rings as 
measured through radial profiles. The best-fit planar ellipses to these points 
are traced by thin, solid lines. 
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FIG. 41. — Rendered points showing the determination of OR positions. 
The NOR and SOR must lie along the upper (opening away from the ob- 
server) and lower (opening toward the observer) medium-grey tubes, respec- 
tively. One possible location of the rings could be the intersection of these 
tubes with the CS hourglass, shown in panels (a)-(b). Panels (c)-(d) show 
the inner edge of the belt (light grey), which is a good fit to the SOR. 



be determined as the intersection of the hourglass with the 
"tube" of all possible z positions for a given point on the ring. 

Figures I41h-b show the results of this exercise. The red 
hourglass from Figure^]is shown in dark grey, while lighter 
dots trace all possible z values for each point along the NOR 
(upper tube) and SOR (lower tube). The hourglass intersects 
both tubes at many points, but fails to do so along the western 
side of the NOR and the entire northern half of the SOR. If we 
instead consider the intersection of each tube with the outer 
edge of the hourglass (using the radial-profile error bars in 
Fig- EH}, the NOR is well fit at all RA.s by points located at 
z' = 1 .37 ± 0.05, while the SOR still fails. 

We consider two possibilities for the SOR. First, were the 
southern lobe oriented with i y ~ 0°, it would cleanly intersect 
the SOR tube. It is only by construction that the hourglass 
has a single axis of symmetry, unfortunately too little of the 
southern lobe was illuminated by echoes to test for such a 
departure. Second, rather than simply lie on or around the 
hourglass, the SOR may result from its intersection with an- 
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FIG. 42. — Probable positions of the ORs (grey) and the planar, elliptical 
fits (black) to each ring, viewed (a) face-on, {b) 60° east, (c) 60° west, and 
(d) orthogonal to the inclination axis. 



TABLE 3 
Best-fit ellipses to the CS rings 



Ring 




y' 




a 


b 







(ly west) 


(ly north) 


(ly) 


(ly) 


(ly) 


(°) 


ER 


0.015 


0.0 


0.0 


0.647 


0.98 


81.1 


NOR 


0.26 


0.04 


-1.36 


1.42 


0.94 


70.5 


SOR 


0.19 


0.06 


1.00 


1.59 


0.92 


-1.1 



other structure. In Figure [38] the dashed line traces a sharp, 
inner edge to the belt, which nearly intersects the southern 
lobe of the hourglass (traced with squares). When the diago- 
nal segment of this boundary is revolved around the hourglass 
axis (light-grey in Fig.l41b-d). it intersects the SOR tube con- 
tinuously at all RA.s at a distance of z' = 1.0 ±0.05 ly from 
the SN. We question whether the SOR is some sort of den- 
sity enhancement resulting from the interaction between the 
hourglass flow and the inner edge of the belt. 

This same belt material makes a poor fit to the NOR, and 
we can not readily identify a similar intersection of observed 
structures. However, echoes from material behind the SN 
could have been missed for a variety of reasons, including 
source confusion and small scattering efficiency. 

The 3-D positions of the ORs are rendered in Figure^] and 
appear reasonably plan ar. W e approximate both rings as pla- 
nar ellipses (Appendix IB. 21 . fit to each distibution of points 
viewed along the CS inclination axis. The resulting parame- 
ters are listed in Table[3] where (x' ,y' ) are the centroid offsets 
from the axis, z' Q is the distance to the SN along the axis, (a, b) 
are the semi-major and minor axes, and <f> is the RA. of the 
major ax is from north. For completeness, the parameters for 
the ER (iy~3 are also listed. 

These fits are indicated in black in Figure |^] and by thin 
white lines in Figure |40| The approximation to the NOR is 
quite good, but the SOR ellipse fails to intersect the points 
that pass just north and west of the ER. This is not suprising, 
since no single ellipse can fully reproduce the observed shape 
of the SOR unless it is non-planar. Since these ellipses are 
consistent with both the CS structure and observed three ring 
nebula, we have used them in the echo renderings throughout 
this work. 
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FIG. 43. — Rendered views of the probable structures containing the observed echoes. Top Row: The CD nebula. Second Row: The NH nebula. Third Row: 
The CS nebula. Fourth Row: The CS hourglass. For these top four rows — Left column: the face-on view of the complete structure; Second column: the western 
half, viewed from the east, sho wing a clear view of the interior; Right columns: face-on and side views of the complete structures in monotone grey, overlaid with 
actual echo points from Figure fl5l For the bottom row — (q) The western halves of the CD and NH, and (r) NH and CS nebulae, viewed from the east, (s) The 
northern halves of the CD and NH, and (f) the NH and CS nebulae, viewed along the inclination axes from the south. See text for colors. 
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9. DENSITY AND MASS OF THE CSE 

9.1. A complete picture of the CSE ofSN 1987 A 

The probable geometry of the structures containing all 
echoes within 25" of SN 1987A are rendered in color in Fig- 
ure ^] and summarized in Table [4] The final column gives 
the approximate total volume of the structure containing each 
echo, found by integrating the radial profiles in Figures [26j 
l32l and[38]along, and revolved 2tt about, the respective incli- 
nation axes. 

The east/west and north/south inclinations are quite consis- 
tent, suggesting all circumstellar material shares a common 
inclination of i x = 40° south and i y = 8° west. The inner two 
structures (CS and NH) both have elliptical cross sections (see 
Figurel43lf-f'). although the ellipticities and orientations differ. 
However, it is unclear whether this particular result is signifi- 
cant. As noted in AppendixIBl \ 2 minimization does not read- 
ily allow for errors in multiple dimensions, yet the measured 
echo positions have finite width and depth. Rather than at- 
tempt much more complicated (and error-prone) techniques, 
such as iteratively reweighted least squares with robust esti- 
mators, we give these results with the understanding that the 
CS and NH ellipticies are most-likely bounded between the 
values reported for the two structures. 

The CS gas is nested neatly within the NH material. The 
extended CS equatorial plane connects to the waist of the 
NH hourglass, forming a continuous structure from 1.5 to 5 
ly from the SN. Tracing NH out from the waist, there are 
inward-pointing spurs above and below this equatorial plane 
that connect to the outer CS wall, which itself connects to the 
belt. The CS belt, equatorial plane, outer wall, and the waist 
of NH are therefore consistent with a single, uniform thick 
waist, extending inward to the CS hourglass. This is shown in 
Figurel43K 

The outer edges of the NH hourglass almost reach the inner 
spur of the CD shell, and indeed the echoes as viewed on the 
sky do connect at early times (Fig. I10i. As shown in Figure 
143b . no material was illuminated along the equatorial plane 
exterior to NH. Unless this r egion w as shadowed from the SN 
by material at smaller radii ( 39.2. U . there is no higher density 
gas outside the waist of NH, suggesting the NH hourglass is 
the pinched waist of the CD. 

The one caveat to this argument is that just above the equa- 
torial plane and outside the NH hourglass are the northwestern 
WFPC2 echoes, whose relationship to the rest of the CD ma- 
terial was unclear in |6] Shown as gold in Figure l43"lw. these 
would have formed a cylindrical band outside the waist of the 
NH hourglass. We do not believe these are image artifacts, nor 
part of a distinct, uniform, and rotationally-symmetric struc- 
ture. Rather, they may indicate local inhomogeneities in the 
undisturbed CSM, or they may simply trace the outermost 
edge of the contact discontinuity. In the following, we will 
exclude them from our analyses. 

Figure |^] shows all circumstellar material, simplified ac- 
cording to the above arguments. The outer structure is peanut- 
like, extended along the poles to r < 28 ly, and narrowly- 
pinched at the waist at r ~ 5 ly. We will henceforth refer 
to this outer shell (traced by the CD and outer NH echoes) as 
the "Peanut." Roughly 35° from the equatorial plane are the 
spurs, extending 3 - 4 ly out from each lobe, perpendicular to 
the Peanut's axis. Unlike the northwestern WFPC2 echoes, 
these spurs do appear at more than one location, but they are 
limited enough in extent that it is unclear whether they are 
uniform features that encircle the entire structure, or are iso- 




Fl G. 4 4. — Re vised cartoon showing the simplified CSE suggested in Figs. 
El E31 and ED Structures are shaded to indicate density, which increases 
with greyscale. Fig. is to scale, and orientation arrows indicate 10 ly. 
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FIG. 45. — Flux density of scattered light per unit volume as a function of 
3-D radius from the SN, measured for each echo in the F612 filter. Points 
are shaded in proportion to their signal-to-noise, which increases with darker 
greyscale. Lines are plotted as follows: j <x r _1 (dashed), j <x r~ 3 (dotted), 
j <x r -4 5 solid. 



lated clumps. 

Inside the equatorial region of the Peanut is a thick (~ 3.5 ly 
in radial thickness, 2 ly in axial length) annular ring (or belt), 
which extends inward until it terminates at the inner CS hour- 
glass. The extent of this waist in z' is unclear. Echoes were 
detected outside the CS hourglass at most equatorial (tt') radii 
between 1.6 and 4 ly, with a vertical (z 1 ) distribution extend- 
ing at least as far as the CS hourglass itself. Still, the bulk of 
material is within 1 ly of the equator. 

9.2. Density of the Circumstellar Material 

To constrain the gas density and dust comp ositio n of the 
CSE, we apply the dust-scattering model from 32.4l using the 
surface brightnesses of each echo discussed in ^JJS] Fluxes 
were only measured in the F470, F612, and F688 filters. Our 
integrated SN spectrum does not extend to 8090A, thus the 
model cannot be applied to the F809 filter without making 
considerable assumptions. We chose not to include measure- 
ments from WFPC2 data for a number of reasons. First, 
echoes were very extended and faint, resulting in large photo- 
metric errors. Second, intercomparing the observations from 
wide-band WFPC2 filters with the narrow-band, continuum- 
only ground-based data is challenging, due to e.g. potential 
line contamination, different filter wavelengths and widths, 
etc. 

Figure l45l shows the volume emissivity j (flux density per 
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TABLE 4 
Summary of Echo Geometries 



Struc- 
ture 


Inclination 3 

ix i y h 

(°) (°) (°) 


Cross Section b 
b/a Axq Ayo 

(ly) dy) 


Shape 


r 

(ly) 


Geometry 

Z' 7T' 

(ly) (ly) 


Att' 

(ly) 


Vtot 

(ly 3 ) 


CS 


41° 


-8° -9° 


0.94 


<o.i 


0. 


Hourglass 


1-2.9 


0.3-2.5 


0.8-1.6 


0.5-1 


20 














Belt 


1.5-2.8 


0.-1.0 


1.5-2.6 


1.5-2.0 


40 














Walls 


1.8-4.2 


0.8-2.7 


1.8-3.6 


2.0-3.0 


120 


NH 


40° 


-7° 12° 


0.82 


0. 


0. 


Hourglass 


3-11 


0.7-7.5 


2.8-7.5 


2-4 


1300 


CD 


40° 


-8° ••• 




<-l 


<-0.1 


Prolate 


10-28 


4-28 


0-16 


4-5 


1.2 x 10 4 



a PYR angles applied in order of roll-pitch-yaw f 85.3.11 . 

b Axo and Ayo are the western and northern offsets of the symmetry axis from the SN, which is at x = 0, y = 0. 

c See Fig. ll2l for coordinate definitions. A-7r' is the approximate width along tt'. Vwi is the inferred total volume of the structure. 



unit volume) of all echoes in F612 as a function of radius. 
Note that 2-D surface brightness was converted to 3-D flux 
density using the line of sight echo depth Az given in Eq. (eq. 
l5ll. T his can be compared with Figure 3 of ICrotts & Kunkell 
(1991). If the density of material is constant, j will decrease 
as r" 1 (Eq. [12] of Sugerman 2003), while a uniform wind 
blowing into a homogenous, low density medium will pro- 
duce an r~ 3 emissivity profile. These are shown in the figure 
as dashed and dotted lines. The outer-NH and CD echoes are 
well-described as constant density, while some NH/CS ma- 
terial may follow an r~ 3 profile. The best fit to the CS and 
inner-NH points follows r" 4 5 , which is suggestive of a dis- 
turbed CSE, su ch as would occur via interacting stellar winds. 

Recall from ^2. 41 that we defined three dust-grain size dis- 
tributions: S, M, and L. We begin with type-L dust. The re- 
sults are shown for pure-carbonaceous (C-only), pure-silicate 
(Si-only), and an LMC-abundance mixture (Ci+S) in Figure 
|46ja). For clarity, densities are plotted logarithmically but po- 
sitions are given linearly, with the horizontal scale changing 
between each vertical line, which roughly delineate CS, NH, 
and CD material. 

The C-only dust yields a close fit between F612 and F688 
for the CS material, but the F470 is uniformly too dense. An 
overdensity implies a large scaling «h, which results from an 
underluminous prediction. Since this dust does not produce 
enough blue flux, the model is too red. The Si-only model 
reconciles the F470 and F688 points, but predicts too much 
flux in F612. Adding carbon to silicates reddens the spec- 
trum, which again pushes F470 densities to be slightly too 
high. Overall, the Si-only is the better fit. 

For the most part, the NH densities vary in a similar fash- 
ion. The C-only model is too red, while the Si-only dust is 
a very good fit. Adding C to Si reddens the model slightly, 
making the F470 densities slightly more discrepant with the 
red points. The Si-only model is again the better fit. Finally, 
the CD echoes are well fit by all three models, with the C+Si 
model having the lowest scatter between filters. 

Figure |46ffo) shows the type-M model results. For the CS 
structures, the C-only model is a reasonable fit to the F470 
and F688 data, yet the F612 flux is too high. Si-only dust 
is overly red, and adding carbon does little to improve the 
fit. Since smaller particles do not scatter enough red flux, the 
CS dust is better fit by a distribution including some large 
grains. The NH material is reasonably-well fit by all three 
models, but the least dispersion between filters occurs when 
silicates are present. Type-M dust is not a good fit for the 
CD, since this model is overluminous in F612 and too blue. 
Although not shown, the fits of the type-S model are poor 



since the spectrum is too blue. 

The CD colors require large grains, and are best-fit by the 
type-L C+Si model. Neither L nor M models are satisfactory 
fits to the NH and CS echoes, which leads us to explore inter- 
mediate grain-size ranges. We tested a number of upper size 
limits between a = 0. 1 -2.0/im, and find that smaller values fit 
better. The CS and inner NH data (r < 5 ly) have the lowest 
dispersion between filters for Si-only, a max = 0.2/im dust. We 
will denote this new size range as ML, as it includes particles 
slightly larger than the type-M dust. Outer-NH data (r > 5 ly) 
are more consistent with Si-dominated dust with larger grains. 
The general trend is that inner echoes are better reproduced by 
smaller, Si-dominated dust, and with increasing distance, the 
g rain sizes and C-content increas es . 

iFisch era. Tuffs, & Volkl ( 120021) have modeled mid-IR emis- 
sion from collisionally-excited dust grains heated in the 
shocks between the H II region interior to the ER and the for- 
ward blast of the SN. They find the emission is best explained 
by small (a < 0.25/im) grains with a Si-Fe or Si-C composi- 
tion. The dust abundance is quite low, which they attribute to 
evaporation from the UV flash and sputtering in the shocked 
gas. This further constrains grain sizes to a < 0.25/im and ex- 
cludes a pure-carbon composition. The dust in both this H II 
region and the CS hourglass was formed from material ex- 
pelled very late in the progenitor's life, thus we expect the dust 
properties to be similar. Indeed, our CS-dust model also fa- 
vors silicate-dominated 11 dust with small grains (a < 0.2/im). 

Our models suggest that dust which formed at later 
times (and is consequently closer to the star) has a smaller 
carbonaceous-component. This can result from a change in 
surface CNO abundances over the star's late-stages of evo- 
lution, since carbon-rich envelopes create carbon-rich dust, 
while oxygen-rich envelopes create silicate-rich dust. Such 
CNO processing is also inferred from early IUE spectra 
(Fransson et al. 1989) of the ER, which show nitrogen and 
oxygen to be overabundant with respect to carbon. 

The average dust density, taken from the best model results 
above is plotted in Figure06jc). The innermost CS hourglass 
material has a relatively-constant gas density of «h = 2 - 3 
cm -3 up to r ~ 1 .6 ly. We note that thi s value is consistent 
with the findings of Burr ows et al.l (119951) . who argued the CS 
hourglass must have a density < 5 cm" 3 to not be observed in 
recombination in WFPC2 images. 

Beyond this position, the material splits into two distribu- 

11 "A stronomica l silicate" dust param eters used in our model 
iDraine & Lee 1984; Weineartner & Draine 2001) are derived from olevine, 
a magnesium-iron-silicate. 
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FIG. 46. — Gas density versus radius for the dust scattering model. Model results are shown for (a) type-L and (b) type-M dust. Each filter is shaded its 
(roughly) corresponding color: F470 in blue, F612 in green, F688 in red. Large circles denote signal-to-noise above 2, and small ones show S/N < 2. The 
absissca changes between columns, each of which contains the structure as noted in the top subpanel. Rows show the results for the dust compositions listed to 
the far left, (c) Average gas den sitie s that match observed echo fluxes to the dust-scattering model. Note that the abscissa is shown in now logarithmic. Point 
colors correspond to those in Fig. 1431 CS material is marked by circles, NH by squares, CD by triangles, and smaller symbols denote measurements with S/N < 2. 
Dotted lines are the best-fits through the upper and lower CS loci, corresponding to n oc r~ 3 (upper) and n oc r^' 5 (lower). Horizontal dashed lines delineate the 
rough density boundaries of the CD gas. 



tions, for which density drops either slowly («h oc r" 3 0±0 - 2 ) or 
rapidly ( r - 4 - 5 ± - 2 ) with radius. Material outside the hourglass 
tends to lie along one of these two distributions, up to a radius 
of roughly 4 ly. That the density varies more steeply than the 
inverse-square expected for a freely-expanding wind suggests 
the mass-loss mechanical luminosity (MvjL,) increased with 
time toward the end of the RSG. 

The densest material traces the waist of the hourglass. The 
hourglass points along the r" 3 profile are those at large tt' for a 
given z', which lie outside the average hourglass position, des- 
ignated by the "extended flux" squares in Figure 1381 Beyond 



r = 2 ly, this profile traces the belt and equatorial plane mate- 
rial. The steeper, r" 4 5 profile includes most hourglass mate- 
rial and extends to include the inner and outer "walls." Allow- 
ing larger dust grains or increasing the carbonaceous content 
of the dust shifts the densities uniformly lower by only 10%, 
thus these results are reasonably robust against the particular 
model chosen. 

The outer CS gas makes a smooth transition to the inner 
NH points around r > 3 ly. The two profiles described above 
also bound the majority of the NH gas up to a radius of 6 ly. 
Excluding a density enhancement between r = 4-5 ly, there 



34 



Sugerman et al. 



is no evidence of a structural distinction between the CS and 
NH gas, which justifies the simplified model in Figure l44l 
Returning to the aforementioned enhancement, there is a thin 
shell or wall at r ~ 4.5 ly, of thickness ~ 1 ly and density «h = 
0.2-0.3 cm" 3 , roughly three times higher than the rest of the 
NH gas. This high-density material is located primarily along 
the waist of the NH hourglass around ir' = 4-5 ly, marking an 
outer edge to the equatorial overdensity. 

The CD density is constant with radius, and bounded by 
0.02 < «h < 0.2 cm" 3 , which are delineated by dashed lines 
Figure I46f c). The greater majority of the NH gas is also 
bounded by these values, confirming the suggestion from Fig- 
ure |^] that the outer NH and CD mark a constant-density 
boundary. Although not a significant trend, the inner CD 
gas (blue conical-shell material), is at slightly higher den- 
sity (n H ~ 0.06) than the red, radial material at larger radii 
(«h ~ 0.04). The green (radial jet) points have the lowest den- 
sity around ~ 0.02 cm" 3 . 

The greyscale shading of structures in Figure l44l reflects 
the differences in density among the many echoing structures, 
with darker grey indicating higher-density material. The pic- 
ture of the CSE will be completed in the following subsec- 
tions, in which we reexamine the Peanut's pinched waist, 
measure the density of dust between the bright echo struc- 
tures, and estimate the total mass of all CSE material. 

9.2.1. The Peanut: a Pinched Waist or Shadowing ? 

According to Figure |^] the CD and NH form a Peanut 
with a narrowly-pinched waist. To interpret this as the con- 
tact discontinuity between the RSG and MS winds, we must 
assume there is no disturbed MS material beyond these struc- 
tures. The alternative is that material is present outside the 
waist of the NH structure, but that it was not illuminated by 
the SN pulse because of significant extinction, or shadowing, 
from equatorial material interior to the Peanut. 

Heuristically, this situation is unlikely. That echoes were 
observed within 0.5 ly of the equatorial plane (Fig. I43K.A 
demonstrates SN light did illuminate this region, so it is dif- 
ficult to argue that the entire equatorial region is shadowed 
from, for example, the ER. An alternative is that all material 
inside NH has shadowed any material outside it. Let us as- 
sume the waist of the Peanut (NH echoes) has a lower density 
than the CS equatorial disk. Extinction scales with density, 
so the waist of the Peanut will cause much less shadowing 
than the material inside it. Otherwise stated, the Peanut itself 
contributes very little to the total optical depth, and therefore 
should not have been illuminated by echoes. Since NH echoes 
were observed at the equator, we conclude that material out- 
side NH should have been illuminated as well, which is not 
the case. 

If, alternatively, the waist of the Peanut has very high den- 
sity and is optically thick, then less light will be scattered to- 
ward the observer than expected in our optically-thin model, 
and we will underpredict the actual dust densities by a large 
margin. We find the NH dust has about the same density as the 
CD shell of the Peanut in the previous section, which means 
the actual NH density must be just the right amount to appear 
to have the same properties as the CD material. This seems 
highly improbable, and leads us to doubt that any strong shad- 
owing has occurred. 

A conservative estimate of the column density of all equa- 
torial dust, by integrating «h x t to a radius of 6 ly in Fig- 
ure 06jc), is Nn ~ 6 x 10 18 cm" 2 . The ER was not ob- 
served in echoes, however an average gas density of 10 4 



cm" 3 (ISo nneborn et alJll997h and observed width of 0.3 ly 
iBurrows*etaniT995t) yields an additio nal column of A^h < 
3 x 10 21 cm" 2 . The total extinction (Sugermar] 120031) at 
6120A for type-ML dust will be A 6 i 20 = 5 x lO" 23 ^ for C+Si 
particles, and 30% smaller for Si-only. The total extinction 
from the SN to NH will be less than 0.2 mags, which will 
only attentuate the SN pulse by 13%. Since the column den- 
sity within the ER is over two orders of magnitude larger than 
in the CSM, any shadowing from the ER would have hidden 
the CS and NH echoes as well. 

If the ER is onl y the ionized s kin of a much more extended 
equatorial disk (Plai t et aT1ll99 5) then the equatorial extinc- 
tion will be correspondingly higher. Only if this disk has con- 
stant density and extends out to 1.5 ly can the equatorial ex- 
tinction approach an optical depth of unity. Again, this seems 
unlikely. The qualitative and quantitative cases for shadow- 
ing are dubious, thus the narrow waist appears to be a genuine 
boundary. 

9.2.2. An Unbiased Measurement of Density 

The previous results are biased, in the sense densities were 
measured only in those regions known to be dense. Any dust 
intersecting an echo parabola will scatter light toward the ob- 
server, thus if any is present between or external to the known 
structures, its signature should be recorded in faint, back- 
ground surface brightness. To probe these regions, we applied 
the dust-scattering model to unbiased photometry, as follows. 

Surface brightness was remeasured in all difference images 
(Fig. I lOt through circular-arc annuli from 0."3 to 50" (or 45" 
until 1990, as explained below) from the SN. Photometry be- 
yond p ~ 25" is intrinsically unreliable, since at these radii we 
made no effort to keep the reference images echo-free. How- 
ever, with the except ion of echo R3 10 (p ~ 45" until 1990), all 
interstellar echoes JXuet al.ll 19941 1 19951) are located beyond 
p = 50", and therefore cause no contamination. The results 
are shown in Figure using the type-L Si+C dust model. 
Note that this figure extends to 300 ly in radius, a factor of 10 
greater than the spatial coverage in Figure l4oT c). 

Panels (a) and (b) differ only in how surface brightnesses 
were measured. Semi-circular annuli, centered at PA. 0° and 
180°, were used in panel (a) to probe uniform, faint variations 
in background flux while smoothing over discrete, non-echo 
signal. Panel (b) shows the results from annuli of arclength 
20°, centered at PA.s 90°, 135°, 180°, 225°, and 270°, which 
better preserve small or non-circular echoes. In both cases, 
the average profile (black points and error bars) is determined 
as the weighted average of densities binned in radius. 

A number of known features are recovered, such as the CS 
material at small radii, the denser waist of the Peanut around 
r = 4 ly, and the bright CD between 1 1 - 22 ly. The average 
density between the CS and CD structures (r = 4 - 22 ly) de- 
creases as r" 1 in panel (a), and r" 2 0±0 5 in panel (b). The exact 
value of the exponent is sensitive to the particular limits in ra- 
dius. Nonetheless, the interior of the CD is filled with diffuse 
material. 

As expected, the semi-circular apertures smoothed over 
most of the CD structure, however, a drop in density around 
r = 30 ly in panel (b) is consistent with its outer boundary. 
External to the CD, the gas density is very uniform with an 
average value of 0.01-0.02 cm" 3 , up to a radius of 80-90 ly, 
at which there is a marginal density peak. These same den- 
sities correspond to that expected from a background-noise- 
only model at r > 30 ly, and are thus only upper limits to the 
actual density in this region. 
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FIG. 47. — Type-L Si+C dust-scattering model densities for all flux within 
50" of the SN, as measured through (a) semi-circular annuli centered at P.A. 
0° and 180°, and (b) 20° annular wedges centered at P.A.s 90°, 135°, 180°, 
225°, and 270°. Signal-to-noise increases with greyscale, and points with 
S/N < 1 are shown as crosses (included since our error estimates are very 
generous). Weighted averages measured through bins of constant A(logr) 
are shown in black. 



9.3. Total Mass of the Circumstellar Environment 

The total mass of gas and dust within the CSE is estimated 
from the average density and volume of each structure. The 
volumes of the probable structures traced out by light echoes 
are measured by integrating the binned radial profiles in Fig- 
ures [33] and [33] Each profile is assumed to represent a 
surface of revolution, and the limits of integration in tt' are 
determined as a convolution of the standard deviation of the 
average profile with the the average spatial extent of each echo 
position, transformed into the cylindrical coordinate system. 
The resulting values are given in the final column of Table |4] 
The volume interior to the CD is estimated a few ways: as 
the volume of the ellipsoid enclosed by the CD, as the volume 
of two cones enclosed by the quasi-linear profile, and as the 
integrated volume under the same profile CFig. I26i. All three 
methods give a consistent result of about 6500 ly 3 . 

The estimated masses of each structure are listed in Table 
|5J yielding total mass in all structures of ~ 1.1 Mq. Since the 
gas-to-dust mass ratio is large (400-600), dust contributes a 
negligible fraction to the total mass. Number density is con- 
verted to mass density by p g = pnunin where [i = 1 .4 is the 
cosmic average value. It is very difficult to quantify all the 
sources of error that enter these values, but we believe the 
tendency will be to underpredict the masses of the structures, 
given conservative measurements of their volume and the un- 
known extent to which dust was not observable, due to un- 
favorable geometry, poor data quality, confusion, or possible 
shadowing. 

10. CONCLUSIONS 

In this work, we have presented the discovery and anal- 
ysis of light echoes illuminating material within 30 ly of 
S N 1987A. These results are also summarized in the first half 
of lSugerman et alJ (|2p05), after which we show that no exist- 
ing formation scenarios can explain this newly-revealed cir- 
cumstellar nebula. 



As noted in the introduction, the outflows from Sk -69° 202 
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Structure 


M/M e 


CS Hourglass 


0.04 


CS Belt 


0.04 


CS Walls 


0.06 


NH Waist 


0.06 


NH Walls 


0.07 


Intra-CD 


0.2 


CD 


1.2 


Total 


1.7 



are one example of a much wider class of bipolar outflows 
seen in many different classes of evolved stars. A small 
sampling of these include the supermassive LBV rj Car, the 
B1.5ab BSG star 25 Sher in NGC 3603, the symbiotic star He 
2-104, the "Etched Hourglass" MyCn 18, and the "butterfly" 
nebula M2-9. The underlying physical processes producing 
these objects are still far from understood. If the shaping of 
bipolar outflows occurs very early in the mass-loss history of 
a star, then the seeds of bipolarity or asymmetry have a strong 
effect on subsequent outflows, but are also erased by those 
same winds. Reconstructing these objects' mass-loss histories 
is very difficult since there are few observational constraints 
extant from their early stages of evolution. Yet such "stellar 
paleontology" is critical to understanding the mechanisms of 
aspherical stellar-mass loss. 

As one researcher recent ly noted at a conferen ce on asym- 
metric planetary nebulae jMeixner et aD|2004), ' 'most as- 
tronomers would gladly give a month's post-doc salary to 
know the three-dimensional structure of a bipolar nebula." 
The light echoes from SN 1987 A have finally provided just 
that, not only for the visible nebula, but additionally for the 
products of stellar mass loss from the last > 3 x 10 5 years of 
the star's life. It is our hope that this work will provide both a 
motivation and a means to readdress first the mechanisms of 
mass loss in Sk -69° 202, and subsequently the phenomena 
of bipolar mass loss in many more systems. 
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APPENDIX 

A. THE DIFIMPHOT DATA PIPELINE 

The difference imaging pipeline requires the following steps: (1) geometric regi stratio n, (2) creating the reference image, (3) 
building the PSF, and (4) PSF-matching and differencing. Step (1) is discussed in 34.3.11 the rest are described below. 

A. 1 . The Reference Image 

A successful implementation of PSF-matched difference imaging will remove all constant sources of flux, leaving the desired 
signal from time-variable sources, and unavoidable contributions from statistical noise, residuals from bright stars, and cosmetic 
defects. Images within a dataset may be differenced against each other, to test for relative changes. Preferably a "reference" (or 
master) image, which represents the median flux of all stars, may be subtracted from all input i mages . In this implementation, 
we combine all images in a given filter to create an "echo-free" reference image, as described in 34.3.21 

A.2. PSF Building 

PSF-matching two images requires high signal-to-noise empirical PSFs with clean profiles well into the wings. Building such 
PSFs can be difficult in sparse and crowded fields. In sparse fields, finding a sufficient number of bright stars to build a PSF 
representing the average stellar profile can be problematic. In crowded fields, isolating bright stars from many nearby, faint 
sources is the dominant challenge. We address both problems by using the PSF-building software from daophot in an iterative 
technique, as follows. 

A candidate list of PSF stars is generated from the reference-star list ( 34.3.2> . using only those members that do not have any 
bright neighbors within > 3 xFWHM pixels. Saturated stars are rejected, but stars with non-linear profiles are kept. The test PSF 
is built from the N brightest linear-profile stars, with non-linear stars used only to supplement the wings. The psf algorithm fits 
an analytic function to each profile and computes the residuals of the star from that function. A weighted average of all profiles, 
normalized to that of the brightest star, is the temporary PSF. This is compared to each input star, and those with residuals that 
deviate from the model are down-weighted in successive iterations. 

The output model PSF is fit to each candidate star in the original list, and up to N of the best fit stars are used as a second input 
list to the PSF-building algorithm. Applying a cut-off to the fit quality is highly effective at removing stars with faint neighbors 
that would otherwise contaminate the wings of the PSF, as well as removing non-linear and blended stars. 

A. 3 . PSF Matching and Differencing 

Both the input and reference images are evaluated for seeing, and the better-seeing (i.e. lower FWHM) of the two is chosen as 
image r. PSF-matching proceeds by first using equation (|8} to build the convolution kernel from the input and reference PSFs. 

Although we take great care to build high-quality PSFs, the noise in the wings can be quite high, since these are the lowest-flux 
components of the radial profile and are therefore the most-poorly determined. As the PSF of a star is roughly gaussian, a fixed 
percentage of the PSF wings is modeled and replaced by the FT of an elliptical gaussian. The particular percentage to replace 
depends on the number and flux of stars used to build the PSF. In most applications, we find no more than 10% of a PSF need be 
replaced, with a well-sampled (~ 20 stars) PSF needing only a 3% replacement. 

Image r is convolved with the kernel k to PS F-match it to the poorer-seeing image i. The photometric scaling c between i and 
r is evaluated from the reference stars ( 34.3.2> present in both images. The difference image is then given by i-c(r*k). It has 
been found that images that are already close in seeing often need to be PSF-matched twice, with the first set of matched images 
as input for a second matching iteration. 

B. NUMERIC MODELING OF DATA 

The generalized problem of finding the parameters a = aj , j = l,M that best-fit the function y(x) = y(x, a) to the N data points 
(yi,Xj) is solved by minimizing the chi-squared merit function 




(Bl) 



where <r, is the weight of y t . In cases in which the function is not linear in x,-, chi-squared minimization requires an iterative 
procedure that relaxes x 2 b y variation of parameters a. Our preferred technique is the use of Levenburg-Marquardt (LM) mini- 
mization JPress et alJl992l) . Briefly, one varies a trial set of parameters a using the method of steepest descent, by computing the 
gradient and Hessian of x 2 ( a )- 

Although very powerful, LM minimization has numerous (yet unavoidable) pitfalls. First, the topography defined by Vx 2 can 
have many local minima into which the fit can relax, missing the absolute minimum. The minimization algorithm is known to 
wander near a minimum if the topography is flat. Complicated functions with many parameters often require that the user hold 
certain ones fixed while varying others, lest the fit reside permanently in local minima. Occasionally, the function can realize 
a minimized \ 2 through an entirely unreasonable set of parameters. Despite these difficulties, the technique is robust and the 
standard of nonlinear fitting. 

Since the datapoints may not have normally-distributed errors, and many fits are to functions of more than one dimension, we 
adopt the bootstrap technique (see iPress et aljfl992l) to estimate the uncertainties in the fitted parameters. Below, we describe 
explicitly a number of non-trivial minimizations used in this work. 
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B . 1 . Gaussians and Moffats 
A moffat is a modified Gaussian, and can be expressed in terms of its peak A, center xq, and FWHM a as 



y = A 



,, 4(2^-1) , 2 ] _/5 
1 + 5 (x-xo) 



(B2) 



where (3 — > oo yields a Gaussian. Any number of moffats can be fit simultaneously, with four parameters per function. 



B.2. Ellipses 

The cartesian expression for an ellipse with center (xo,yo) and semi-major and minor axes (a,b) is not a form condusive to 
LM fitting for m any reasons, not least of which is that it is particularly complicated to account for rotations by some angle a 
(see Figure lB48h ). Rather, the ellipse can be expressed in parametric form as (x = acos£, y = fosin£) with parameter £. To fit 
an ellipse to a locus of points, we minimize the geometric distance between the ellipse and the data as follows. Data are shifted 
by (— Xq,- yo) and rotated by -a, then for each datum, the value of £ is found that minimizes the point's distance from the ellipse 
defined by (a,b). This specific minimization can be accomplished by simple numeric bracketing, and as £ is only used to find 
points on the ellipse, not to define it, it is not fit by the LM algorithm. The \ 2 merit function is the sum of the distance of all 
transformed data from the ellipse, and is minimized by varying the five parameters {(xo,yo),(a,b),a}. 

Explicitly, this is a two-dimensional fit for the distance function d = c/(x, ,a) = yJXf + Yf, where 

Xi = xi - (xq + a cos a cos + b sin a sin £,) 
Yi = yi - (yo - a sin a cos £; + b cos a sin £, ) . 

The absolute minimum dist ance of all data to the ellipse is zero, thus we wish to find the parameters such that c/(x, ,a) « 0. In the 
relevant analog to eauation lBll \ 2 = S;^(X;,a) 2 . N.B. Here, % 2 is a measure of the total distance of all points from the ellipse, 
and does not carry statistical interpretations such as absolute goodness-of-fit. Rather, -\/x 2 /(N— 1), where N is the number of 
points to be fit, is the root-mean-squared (RMS) scatter of the fit. 

B.3. Three-Dimensional Lines 

Finding the orientation and position of a 3-D line that fits data is equivalent to finding the transformation that aligns the data 
with a fixed line. Define that transformation as the vector offset (xojJOjZo) an< ^ pitch-yaw-roll (PYR) rotation through angles 
(9, <p, i[>)- As for the ellipse, we seek to find the six parameters that minimize the total distance from all transformed points to the 
Z axis. The same caveat discussed above, regarding the statistical interpretation of \ 2 , applies. 

B.4. Cones and Frusta 

As with the ellipse, the traditional cartesian definition of a cone is not condusive to LM minimization. Rather, we implement 
a novel combination of the previous two subsections. A cone can be considered as a series of nested ellipses aligned along a 
given vector n, where the axes of each ellipse are functions of the distance of its center from some point xo. If the ellipses have 
non-zero major and minor axes at xp , the figure is a truncated cone or frustum. The distance from the apex to the truncated top 
edge is defined as h (see Figure lB48b ). 

Proceeding as before, the data are rotated and translated to fit them to a cone aligned along the z axis. The major and minor 
axes are defined to be a(z) = ao(\z\ + h), b(z) = bo(\z\ + h), where h is defined above. Defined in terms of |z|, an hourglass (two 
cones reflected about the x—y plane) can al so be fit. 

Consider the point P as shown in Figure lB48t >. In truth, it is a distance d from the cone, however calculation of the point R, 
or its projection C on the z axis, is a complicated function with extremely complicated partial derivates. Rather, note that the 
triangles AAQB and APRQ are similar, and if the cone has a fixed half-opening angle (ZBAQ = a), d/l = cos(ZQPR) is a fixed 
ratio, thus minmizing / is equiv alent to minimizing d. The application of this is simple, since the distance / can be measured 
using the technique from SIB. 21 Note however that the half-opening angle a is fixed only for a cone with circular cross section 




(ao = bo). For the general case of an elliptic cross section, d/l is not constant with angle about the z axis. However, restricting 
this technique to small eccentricities (a < b), d/l is approximately constant, and this technique is still valid. 

Requiring that the axis of the cone pass through the origin, the center Xo = kn, where h is the unit vector k rotated by the PYR 
angles (ip z ,8 x ,(j)y). This cone (or frustum) can be fit by seven parameters: three rotation angles, two ellipse axes, h and k. For 
a general center xo, the fit requires nine parameters. However, in the adopted PYR rotation convention, ty z has no effect for a 
circular cross section, for which case the above fits are reduced to five and seven parameters, respectively. 

C. VISUALIZATION TRANSFORMATIONS 

A very large number of projectio ns exi st to represent a 3-D object on the two-dimensional image plane of a piece of paper. 
Three of these are shown in Figure IC491 A clinographic projection views an object from an angle <f> with each point on the 
cube projected onto the image plane by a parallel line. The cube is now seen including the top face, and the size h of the cube 
is preserved. An orthographic projection first rotates the cube by the angle <f> and then projects the image onto the plane with 
orthogonal lines, such that the side h now has a length hcoscf>. Such "parallel projections" are difficult to interpret visually since 
we do not observe 3-D structures this way. Rather, our eyes apply conical or perspective projections, by which the relative sizes of 
points diminish with increasing distance from the observer. Since parallel lines project to an imaginary apex, zonal relationships 
appear slightly distorted, but are also easier to understand. 

To view a 3-D object, one first applies orthogonal transformations (positional shifts, rotations, scalings), then the resulting 
object is projected onto the image plane using either a parallel or perspective transformation. The general case provides for: a 
linear shift in three dimensions, rotations about three axes, a uniform scaling, followed by a dimetric and perspective projection. 
Consider the image plane as a vertical pane of glass, with the object held in front of the pane by a metal rod, connected orthogo- 
nally to the glass. The dimetric projection is an orthographic projection applied in two dimensions, and establishes where along 
the glass the observer is located with repsect to that rod. The perspective projection corrects for the distance of the observer 
behind the glass. 

All of these can be easily specified as 4 x 4 matrices. A rotation R, shift S, and translation T are given by 

1 ON 
10 
10 

K Ax Ay Az 1, 

where for R we have left the sin and cos coefficients out since they vary with the particular convention used (i.e. Euler angles 
versus pitch-yaw-roll). A dimetric projection Pj extrudes projected surfaces in the direction of a viewing vector v. A perspective 
projection P p is defined by the parameter a, where a — > oo reduces to an orthographic projection with <f> = 0. They are given by 

P d = 

To apply a perspective projection, the position vector (x,y,z, 1) is multiplied by P^ and then normalized by the fourth component, 
i.e. a perspective-transformed point is given by 

ax ay z 
z+a 1 z+a' z+a 

Any combination of orthogonal and perspective transformations can be applied by finding the matrix product of all relative 
transformations, and then applying that new matrix to each 3-D datum, loaded into a 4-vector as above. The projected position 
onto the image plane is given by the first two components of the final vector, the third component fills the z-buffer to store depth 
information, and the fourth component returns to 1 . 





( a 


b 


c 0\ 




As 






R = 


d 

8 


e 
h 


/ o 

i 


s = 


s 
s 


;) 


-( 




VP 





\) 




\0 







i 














1 








'xl V z 


-Vyh'z 





















The 3-D Map of SN 1987A 



39 



In our implementation of the 3-D rendering described in 35.3.21 we find the dimetric projection is never needed. Most render- 
ings in this work have been generated with a = 10, which creates a shallow perspective that aids the eye in understanding the 
correct orientation of the coordinate box, but has a small effect on zonal relationships. 
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